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1.  INTRODUCTION 


Composite  materials  have  been  increasingly  used  during  the  last  decades,  in 
order  to  lighten  structures  in  fields  like  aeronautics  and  space.  Two  steps  are 
essential  in  the  process  of  taking  advantage  of  this  kind  of  materials:  design 
and  optimization. 

Optimization  of  composite  structures  is  a  recent  issue,  because  both 
optimization  techniques  and  composite  structures  have  been  developed  during 
the  last  few  decades  and  therefore,  the  conjunction  of  both  of  them  is  even 
more  recent.  Composite  materials  being  an  expensive  but  efficient  technology 
to  get  minimum  weight  structures,  it  is  logical  to  make  an  attempt  to  find 
out  how  to  design  properly  optimum  laminated  composite  plates  with  no 
reduction  in  their  strength.  A  general  scheme  is  depicted  in  Figure  1. 

Since  many  kinds  of  ground  and  air  vehicles  have  rectangular  plates  as  a 
common  structural  element,  an  increasing  demand  for  improved  structural 
efficiency  in  such  applications  has  resulted.  Composite  materials  offer  a 
number  of  advantages  other  than  their  high  stiffness  to  density  values:  for 
example,  their  capability  of  being  tailored  by  orientation  of  the  filaments  in 
the  various  layers  in  order  to  optimize  the  desired  structural  behavior. 

Actually,  although  there  is  a  lot  of  literature  related  to  plate  analysis,  the 
most  modem  references  lack  adequate  information  which  could  allow  a 
designer  to  tailor  or  synthesize  an  optimal  design.  The  reason  for  that  can  be 
found  in  two  aspects:  first,  the  difficulty  of  accurately  assessing  the  level  of 
strains  and  stresses  in  any  point  of  a  laminated  composite  plate  and,  secondly, 
the  need  of  implementing  an  optimization  procedure  able  to  find  the 
minimum  weight  structure. 

Fortunately,  the  current  development  of  numerical  techniques  and  the 
existence  of  powerful  computers  make  possible  a  solution  to  the  two  problems 
mentioned  above. 
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The  interest  of  the  optimization  of  laminated  composite  plates  is  focussed 
on  three  different  fields:  in-plane  loads,  transverse  loads,  and  buckling.  It  is 
noticeable  by  reviewing  the  literature,  that  the  problem  dealing  with  in-plane 
loads  has  already  been  studied,  and  a  number  of  technical  publications  and 
software  packages  can  be  found. 

In  the  present  report,  an  optimization  of  rectangular  laminated  composite 
plates  subjected  to  transverse  loads  and  buckling  is  carried  out.  Plates  (both 
thin  and  thick)  and  sandwich  panels,  subjected  to  uniformly  distributed  and 
point  loads,  are  analyzed  in  the  chapter  about  transverse  loads,  considering 
clamped  and  simply  supported  as  boundary  conditions.  Uniform  uni-,  bi-axial 
compression,  shear,  and  combined  loads  are  also  studied  in  the  section 
devoted  to  buckling. 


1.1  LITERATURE  REVIEW 


Optimization  procedures  applied  to  laminated  composite  plates 


The  earliest  attempt  in  composite  optimization  seems  to  be  due  to  Foye  [1], 
who  studied  the  minimum  weight  optimum  design  of  laminated  for  strength 
and  membrane  stiffness.  Multiple  in-plane  loading  conditions  were 
considered,  and  a  random  search  method  was  used  to  find  ply  orientation 
angles,  such  that  the  strength  and  stiffness  requirement  would  be  satisfied 
with  the  smallest  number  of  plies.  Another  procedure  for  the  optimum  design 
of  laminates  was  reported  by  Waddoups  [2].  Minimum  weight  designs  are 
obtained  by  considering  strength  constraints  under  multiple  distinct  loading 
conditions.  Either  the  Tsai-Hill  or  the  maximum  strain  criterion  may  be  used, 
and  all  laminae  are  assumed  to  behave  linearly  to  failure. 

Waddoups  and  co-workers  [3]  reported  a  structural  synthesis  capability  for 
a  class  of  anisotropic  plate  structures.  The  basic  concepts  of  structural 
synthesis  were  developed  by  Schmit  and  his  associates  [4,5].  This  technique 
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was  successfully  applied  to  the  optimum  design  of  integrally  stiffened 
cylindrical  shells  by  Kicher  [6]  and  Morrow-Kicher  [7].  The  most  general 
problem  dealing  with  this  method  was  formulated  by  Waddoups  and  his 
workers  [3]  and  involves  21  design  variables,  45  distinct  failure  modes,  and  3 
independent  loading  conditions.  Kicher  and  Chao  [8]  reported  the 
development  of  a  structural  optimization  capability  for  stiffened  fiber- 
composite  cylinders.  Multiple  load  conditions  were  considered  and  strength, 
as  well  as  buckling  failure  modes,  were  guarded  against.  The  design  variables 
include  stiffener  dimensions  and  spacing,  fiber  volume  content,  and  ply 
orientation  angles.  In  these  last  two  references,  [3]  and  [8],  the  optimization 
problem  is  cast  in  mathematical  programming  form  and  then  a  SUMT-type 
method,  based  on  the  Fiacco-McCormick  interior  penalty  function  and  the 
Davidson-Fletcher-Powell  algorithm  for  unconstrained  minimization,  is 
employed  to  obtain  numerical  results.  Verette  [9]  extended  the  laminate 
optimization  procedure  reported  by  Waddoups  [2]. 

Morz  [10]  tried  to  obtain  the  optimal  reinforcement  using  strength  as  a 
design  criterion.  Bryzgalin  [11]  and  Love  and  Melchers  [12]  used  a  stiffness 
criterion  to  get  an  optimal  design  of  constant  thickness  composites.  Lai  and 
Abenbach  [13]  used  a  direct  search  procedure  to  obtain  an  optional  design  for 
minimum  tensile  stress  at  the  interface  in  a  layered  structure  subjected  to  time 
harmonic  and  transient  loads. 

Khot  [14]  suggested  an  efficient  optimization  technique  based  on  strain 
energy  distribution  and  a  numerical  search  for  the  minimum  weight  design  of 
structures.  The  procedure  takes  into  account  multiple  loading  conditions  and 
displacement  constraints  on  the  structure.  Schmit  and  Farshi  [15,16]  presented 
a  method  for  minimum  weight  design  of  symmetric  composite  laminates 
subjected  to  multiple  in-plane  loading  conditions.  Hayashi  [17]  optimized  the 
fiber  volume  fractions  for  columns  and  orientation  angle  for  plates  and 
cylinders  for  corresponding  buckling  strength. 

Chao  and  others  [18]  determined  the  optimum  fiber  orientation  for  a 
symmetric  orthotropic  composite  laminate  with  in-plane  loading.  The  fiber 
orientation  of  each  ply  is  treated  as  a  design  variable.  Hayashi  [19]  and  Bert 
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[20]  performed  a  similar  analysis.  Housner  and  Stein  [21]  calculated  the 
optimum  fiber  direction  of  graphite-epoxy  sandwich  panels  under  axial 
compression,  assuming  that  the  angle  of  all  the  plies  was  the  same.  Hirano 
[22]  optimized  the  buckling  load  of  laminated  plates  under  uniaxial  and 
biaxial  compression. 

Bert  [23]  presented  a  technique  for  obtaining  the  optimal  laminate  design 
for  a  thin  plate  which  consisted  of  multiple  layers  of  equal  thickness 
composite  laminae,  the  design  criterion  being  maximization  of  fundamental 
frequency.  By  means  of  this  technique,  Bert  [24]  analyzed  clamped 
composite-materials  plates,  Adali  [25]  reported  the  design  of  shear- 
deformable  antisymmetric  angle-ply  laminates,  and  Pedersen  [26,27]  studied 
the  sensitivity  analysis  and  optimal  design  for  laminates. 

Joshi  and  Iyengar  [28]  studied  the  minimum  weight  design  of  composite 
plates  subjected  to  in-plane  and  transverse  loads  treating  fiber  orientation  and 
thickness  of  plies  as  design  variables.  Soni  and  Iyengar  [29]  reported  their 
studies  on  optimum  design  of  damped  laminated  plates.  In  a  design  paper 
Joshi  and  Iyengar  [30]  presented  their  studies  on  optimum  design  under  in¬ 
plane  loads  with  multiple  design  variables. 

Khot  [31]  presented  a  computer  program  (OPTCOMP)  for  analyzing 
composite  structures  subjected  to  in-plane  and  transverse  loads.  The 
recurrence  relations  used  are  remarkable,  but  this  technique  does  not  account 
for  interlaminar  stresses.  McKeown  [32],  Starnes  and  Haftka  [33],  and  Schmit 
and  Mehrinfar  [34]  reported  studies  on  structural  synthesis  systems  dealing 
with  optimization  of  composite  structures,  but  the  design  flexibility  is  limited. 
Sobieszczanski-Sobiesky  [35]  described  a  design  system,  but  he  did  not  fully 
consider  the  interaction  of  all  the  design  variables. 

Park  [36]  analyzed  the  optimal  design  of  simple  symmetric  laminates 
under  the  first  ply  failure  criterion.  Donaldson  [37]  reported  simplified 
weight-saving  techniques  for  composite  panels.  Wurzel  [38]  studied  optimal 
design  of  bidirectional  composites. 
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Massard  [39]  reported  a  computer  sizing  study  dealing  with  optimization  of 
composite  laminates  subjected  to  in-plane  loads.  Maksimovic  [40]  developed  a 
procedure  for  obtaining  near-minimum  weight  design  of  large  composite 
structures  subjected  to  statically  applied  loads  by  means  of  the  method  of 
inscribed  hyperspheres.  Tsai  [41]  described  a  computer  program 
(LAMRANK)  to  optimize  and  rank  composite  laminates  subjected  to  in-plane 
loads.  Watkins  and  Morris  [42]  presented  a  multicriteria  objective  function 
optimization  scheme,  with  application  to  laminated  plates  subjected  to  in-plane 
and  transverse  loads. 

Iyengar  [43]  obtained  an  optimum  weight  design  of  hybrids  laminates  of 
fiber-reinforced  composites  made  of  boron/epoxy  and  aramid/epoxy. 
Pedersen  [44]  studied  how  an  orthotropic  material  must  be  oriented  in  order 
to  obtain  extreme  energy  density. 

A  thorough  study  of  buckling  problems  for  laminated  composite  plates 
was  presented  by  Leissa  [45].  Muc  [46]  obtained  three  types  of  optimal  angle 
in  closed  analytical  form  for  laminated  rectangular  plates  under  uniaxial  and 
biaxial  compression.  These  directions  are  functions  of  the  material  properties, 
the  geometrical  properties  and  the  type  of  buckling  mode. 

After  reviewing  the  literature  related  to  optimization  of  composite  plates,  it 
is  noticeable  that  laminated  composite  plates  under  in-plane  loads  can  be 
optimized  by  following  the  papers  described  above.  There  is  no  difficulty  in 
any  case,  and  comercial  software,  from  personal  computers  to  macro¬ 
systems,  is  available. 

In  the  chapter  about  laminated  composite  plates  subjected  to  transverse 
loads,  none  of  the  papers  mentioned  above  considered  interlaminar  stresses, 
and  there  are  no  solutions  available  in  the  literature  for  different  aspect 
ratios,  boundary  conditions,  and  types  of  loading.  Owing  to  the  complexity  of 
the  problem,  laminated  composite  plates  subjected  to  transverse  loads  were 
analyzed  by  means  of  the  Reissner-Stavsky  theory  [47]  and  therefore, 
interlaminar  stresses  were  neglected. 
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Though  thousands  of  papers  dealing  with  buckling  of  plates  have  been 
written  up  to  now,  there  is  no  publication  related  to  the  optima  solutions  to 
different  problems  dealing  with  one-dimensional  and  two-dimensional 
laminated  composite  plates. 


The  finite  element  method  applied  to  laminated  composite  plates 


In  the  process  of  optimizing  a  composite  structure,  not  only  an 
optimization  procedure  is  needed,  but  a  finite  element  formulation  is  also 
required.  The  following  paragraphs  describe  a  brief  review  about  the  finite 
element  formulation  applied  to  laminated  composite  plates. 

The  classical  lamination  theory  based  on  the  Kirchoff  [48]  hypothesis  was 
the  first  attempt  in  the  development  of  laminated  plate  theory.  This  theory  is 
not  computationally  efficient  from  the  simple  finite  element  formulation  point 
of  view  [49]  neither  is  it  suitable  for  composite  structures  analysis,  owing  to 
its  simplifying  assumptions,  the  most  important  of  which  are  the  neglect  of 
the  transverse  shear  deformations,  and  the  transverse  normal  stresses. 

Reissner  [50]  and  Mindlin  [51]  were  the  first  to  provide  first-order  shear 
deformable  theories  based  on  the  thin  plate  assumptions  for  variation  of 
stresses  and  displacements  through  the  thickness  of  the  plate,  respectively. 
Studies  of  low  order  elements  i.e.,  3-noded  triangles,  4-noded  and  8-noded 
quadrilaterals  showed  violent  stress  oscillations.  Some  techniques,  such  as 
reduced  and  selective  integrations  [52-55],  modified  shear  strain  methods 
[56],  and  hybrid  and  mixed  methods  [57],  succeeded  in  generating  efficient 
elements.  But,  even  so,  these  numerical  analysis  have  certain  limitations:  the 
transverse  shearing  strains  are  assumed  constant  through  the  plate  thickness, 
and  a  fictitious  shear  correction  coefficient  is  introduced. 

Lo,  Christensen  and  Wu  [58-59],  and  Reissner  [60]  presented  a  theory  for 
plates  based  on  an  assumed  higher  order  displacement  field.  Kant,  Owen,  and 
Zienkiewicz  [61]  presented,  for  the  first  time,  a  C°  plate  bending  finite 
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element  formulation  of  a  higher  order  theory.  This  element  is  a  9-noded 
Lagrange  quadrilateral  and  has  six  degrees  of  freedom  per  node.  Reddy  [62] 
presented  a  higher  order  shear  theory  in  which  in-plane  displacements  are 
expanded  as  cubic  functions  of  the  thickness  coordinate,  while  the  transverse 
deflection  is  kept  only  a  function  of  x  and  y. 

Recently,  Pandya  and  Kant  [63]  presented  a  refined  higher  order  C°  plate 
bending  element  in  which  the  transverse  deflection  is  expressed  as  a  quadratic 
function  of  the  thickness.  This  element  is  a  9-noded  Lagrange  quadrilateral. 
Lakshminarayana  and  Ramani  [64]  studied  some  shear  flexible  triangular 
laminated  composite  plate  finite  elements.  Prathap  and  Somashekar  [65] 
studied  a  8-noded  laminated  anisotropic  finite  element.  Tessler  [66]  analyzed  a 
higher  order  theory  in  which  the  transverse  deflection  is  expressed  as  a 
quadratic  function  of  the  thickness.  Miravete  [67]  presented  a  finite  element 
formulation  applied  to  three-dimensional  laminated  composite  plates,  based 
on  [62]  and  the  penalty  function  method  [68]. 

A  more  detailed  description  of  the  finite  element  formulations  mentioned 
above  can  be  found  in  [69]. 

In  this  report  some  finite  elements  have  been  used,  the  formulation  varying 
in  function  on  the  type  of  problem  to  solve. 


1 .2  PROJECT  SUMMARY 


In  this  work,  an  attempt  is  made  to  optimize  rectangular  laminated  composite 
plates.  Laminate  thickness  and  fiber  orientation  are  treated  as  design 
variables.  Since  optimization  of  composite  plates  subjected  to  in-plane  loads 
has  already  been  reported,  the  present  study  is  focussed  on  transverse  loads 
and  buckling. 

In  Section  2  the  symbols  and  the  notation  used  in  this  report  are  described. 
Basically,  this  reports  follows  the  notation  introduced  by  Tsai  and  Hahn  [71] 
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in  those  sections  related  to  design  of  composites,  the  symbology  and  notation 
used  by  Khot  [31J  in  those  chapters  dealing  with  optimization,  and  finally,  the 
notation  introduced  by  Zienkiewicz  [68],  in  the  paragraphs  devoted  to  the 
finite  element  method. 

Optimization  of  laminated  composite  plates  subjected  to  transverse  loads  is 
treated  in  Section  3.  First,  the  equations  of  the  analysis  are  formulated.  Since 
strength  has  been  the  design  criterion,  a  failure  criterion  must  be  chosen. 
Next,  recurrence  relations  dealing  with  the  optimization  procedure  are 
described. 

Once  the  analytical  formulation  has  been  reported,  the  variable  thickness 
problem  is  analyzed  by  means  of  a  2-D  model.  Through-thickness  stress 
distributions  are  determined  for  different  cases,  and  the  influence  of 
geometric  parameters  and  the  applied  load  on  the  strength  of  the  plate  is  also 
analyzed.  The  failure  mechanisms  reported  on  experimental  testings  are 
described  next.  Finally,  verification  of  the  theoretical  model  by  means  of  a 
number  of  experimental  testings  is  carried  out. 

The  final  part  of  the  chapter,  related  to  optimization  of  laminated 
composite  plates  subjected  to  transverse  loads,  is  devoted  to  describe  the 
numerical  results.  The  6-component  stress  tensor  is  accounted  for.  Normal 
and  shear  interlaminar  effects  are  present  in  the  analysis  and,  therefore, 
elastic  constants  and  strengths  must  be  determined.  Optima  configurations  of 
variable-thicknesses  and  angles  of  orientation  of  fibers  are  presented  for  one- 
and  two-dimensional  plates.  Weight  savings  are  also  reported. 

Section  4  deals  with  the  optimization  of  laminated  composite  plates 
subjected  to  buckling.  First,  the  analytical  formulation  is  described. 
Numerical  results  are  presented  for  uniform  uniaxial  compressive  loads  (one¬ 
dimensional  plates)  and  uniform  uniaxial,  biaxial  compressive,  shear,  and 
combined  loads  (two-dimensional  plates).  The  optimization  criterion  used  has 
been  the  critical  compressive  force. 

Concluding  remarks  are  outlined  in  Section  5. 
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2.  SYMBOLS  AND  NOTATION 

a  Length  of  the  plate;  or 

nodal  displacement  vector. 

a '  Relative  nodal  displacement. 

ae  Nodal  displacement  vector  of  the  finite  element  e. 

B  Geometric  matrix  used  in  the  finite  element  method;  B  =  LN. 

b  Width  of  the  plate. 

Qj  Stiffness  matrix  in  the  generalized  Hooke’s  law;  in  concentrated 

notation,  i,j=l,2,3,4,5,6. 

c  Strain-displacement  matrix,  or 

core  thickness  in  a  sandwich  panel. 

D  Matrix  of  elastic  constants. 

ey  Strain  energy  of  the  jth  layer  of  the  element  i. 

ejj  Strain  energy  density  of  the  jth  layer  of  the  element  i. 

e'jj  Strain  energy  density  of  the  jth  layer  of  the  element  i. 

e'ij  Relative  strain  energy  of  the  jth  layer  of  the  element  i. 

Fjj,F,  Strength  parameters  in  stress  formulation  of  the  quadratic  failure 
criterion. 

fe  Nodal  forces  in  the  finite  element  e. 

K  Stiffness  matrix  of  the  whole  plate. 


9 


K' 


Relative  stiffness  matrix  of  the  whole  plate. 

Kb  Stiffness  matrix  of  the  whole  plate  due  to  bending  stresses. 

Ks  Geometric  matrix  of  the  whole  plate. 

ke  Elemental  stiffness  matrix  . 

1  Diagonal  of  a  plate, 

li  Surface  area  of  the  element  i. 

M  Applied  bending  moment. 

N  Shape  function  matrix. 

Nx  Applied  compressive  force  in  x-direction  per  unit  width. 

Nx  Cr  Critical  compressive  force  in  x-direction  per  unit  width. 

Ny  Applied  compressive  force  in  y-direction  per  unit  width. 

Ny  cr  Critical  compressive  force  in  y-direction  per  unit  width. 

NXy  Applied  compressive  force  in  xy-  or  12-plane  per  unit  width. 

NXy  cr  Critical  compressive  force  in  xy-  or  12-plane  per  unit  width. 

P  Applied  transverse  load. 

Q  Reduced  stiffness  matrix  for  plane  stress. 

qe  Distributed  forces  in  the  finite  element  e. 

Qyy  Term  of  the  reduced  matrix  for  plane  stress. 

Qyy  =  Ey  /  (l-vi2v2i). 

r  Strength/stress  ratio  or  strength  ratio. 

S  Positive  pure  shear  strength  in  the  xy-  or  12-plane  of  a  ply. 


10 


S'  Negative  pure  shear  strength  in  the  xy-  or  12-plane  of  a  ply. 

t  Total  laminate  thickness. 

ti  Difference  between  the  total  laminate  thickness(t)  and  the  minor 

thickness  of  a  tapered  plate. 

lij  Thickness  of  the  jth  layer  of  the  ith  element, 

u  Deflection  in  the  x-direction. 

V  Strain  energy  of  a  single  plate  element  subjected  to  in-plane 
stresses. 

Vb  Strain  energy  of  a  single  plate  element  due  to  plate  bending. 

Ve  Volume  of  the  finite  element  e. 

Vij  Volume  of  the  jth  layer  of  the  ith  element. 

Vs  Strain  energy  of  a  single  plate  element  due  to  in-plane  stresses, 

v  Deflection  in  the  y-direction. 

w  Deflection  in  the  z-direction. 

X  Uniaxial  tensile  strength  of  a  ply  along  the  x-axis. 

X'  Uniaxial  compressive  strength  of  a  ply  along  the  x-axis. 

Y  Uniaxial  tensile  strength  of  a  ply  along  the  y-axis. 

Y'  Uniaxial  compressive  strength  of  a  ply  along  the  y-axis. 

Z  Uniaxial  tensile  strength  of  a  ply  along  the  z-axis. 

Z'  Uniaxial  compressive  strength  of  a  ply  along  the  z-axis. 

aij  Thickness  of  the  jth  layer  of  the  ith  element  normalized  to 

maximum  thickness. 
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Strain  generated  by  a  virtual  displacement. 

Displacement  generated  by  a  virtual  displacement. 

Vector  of  nodal  deflections. 

Virtual  displacements  in  the  finite  element  e. 

Strain  component. 

Relative  strain  component. 

Initial  strains. 

Angle  between  x-axis  and  fiber  or  principal  axis  of  layer. 
Scaling  parameter. 

Lagrangian  parameter;  or 

coefficient  used  to  formulate  a  eigenvalue  problem. 
Number  of  cycle  of  iteration. 

Angle  of  variation  of  thickness. 

Mass  density  of  the  jth  layer  of  the  ith  element. 

Stress  component. 

Relative  stress  component. 

Initial  stresses. 


One-dimensional  thin  plate. 
One-dimensional  thick  plate. 

One-dimensional  sandwich  plate. 
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One-dimensional  simply  supported  plate. 


One-dimensional  clamped  plate. 
One-dimensional  cantilever  plate. 

Two-dimensional  thin  plate. 

Two-dimensional  thick  plate. 

Two-dimensional  sandwich  plate. 
Two-dimensional  simply  supported  plate. 

Two-dimensional  clamped  plate. 

One-dimensional  plate  subjected  to  a 
uniform  transverse  load. 

One-dimensional  plate  subjected  to  a  point 
transverse  load. 

Two-dimensional  plate  subjected  to  a  uniform 
transverse  load. 


Two-dimensional  plate  subjected  to  a  point 
transverse  load. 
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One-dimensional  simply  supported  plate. 


One-dimensional  clamped  plate. 


One-dimensional  cantilever  plate. 


Two-dimensional  simply  supported  plate. 


Two-dimensional  simply  supported  plate. 


Two-dimensional  plate  subjected  to  a  uniaxial 
uniform  compression  load. 


Two-dimensional  plate  subjected  to  a  biaxial 
uniform  compression  load  (Ny=Nx/2). 


14 


nraim 


Two-dimensional  plate  subjected  to  a  biaxial 
uniform  compression  load  (Ny=Nx). 


Two-dimensional  plate  subjected  to  a  biaxial 
uniform  compression  load  (Ny=Nx*2). 


Two-dimensional  plate  subjected  to  a  uniform 
in-plane  shear  load . 


Two-dimensional  plate  subjected  to  a  uniaxial 
uniform  compression  and  in-plane  shear  load  . 
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3.  OPTIMIZATION  OF  LAMINATED  COMPOSITE  PLATES 
SUBJECTED  TO  TRANSVERSE  LOADS 


This  section  deals  with  transverse  loads.  The  study  of  the  optimization  of 
laminated  composite  plates  subjected  to  this  type  of  load  is  extremely 
complex. 

On  the  one  hand,  the  whole  stress  tensor  must  be  considered:  the  three  in¬ 
plane  stress  components  are  present  owing  to  the  bending  effects,  and  so  are 
the  two  interlaminar  shear  stress  components,  owing  to  the  shear  effects. 
Finally,  the  interlaminar  normal  stress  component  must  also  be  considered 
because  of  the  variable  thickness,  as  equilibrium  equations  predict. 

On  the  other  hand,  the  minimum  weight  plate  must  be  found.  Therefore, 
accuracy  and  efficiency  should  be  combined  in  order  to  analyze  the  structure 
properly  and  then,  to  solve  the  optimization  problem  by  means  of  an  iterative 
procedure. 

Thus,  first  of  all,  the  analytical  formulation  is  exposed.  General  procedures 
for  analyzing  the  structure,  predicting  the  laminate  failure,  and  optimizing 
the  laminate  composite  plate  are  described. 

Secondly,  the  variable  thickness  problem  is  analyzed.  The  behavior  of  a 
bidimensional  plate  with  variable  thickness  is  extremely  difficult  to 
understand,  because  it  is  a  3-D  problem,  and  because  bending,  shear,  and 
variable  thickness  effects  appear  simultaneously.  Therefore,  a  one¬ 
dimensional  plate  with  variable  thickness  is  studied  here  by  means  of  a  2-D 
plane  strain  finite  element  model. 

Finally,  once  the  failure  mechanisms  are  identified,  and  the  2-D  plane 
strain  model  has  been  verified  by  means  of  an  experimental  study,  a  number 
of  results  are  presented. 
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One-dimensional  laminated  composite  plates  are  studied  by  applying  the  2- 
D  plane  strain  model  mentioned  above.  Optima  configurations  and  weight 
savings  are  reported  for  different  kinds  of  plates,  types  of  loading,  and 
boundary  conditions. 

Two-dimensional  laminated  composite  plates  are  also  analyzed.  Since  the  2- 
D  plane  strain  model  can  only  be  used  for  the  one-dimensional  case,  a  more 
general  model  based  on  a  higher  order  shear  theory  is  used.  This  method  is 
efficient  and  accurate,  not  only  because  the  analysis  is  carried  out  very  fast, 
but  also  because  the  required  number  of  nodes  of  the  mesh  is  not  large,  and 
the  results  from  the  shear  deformation  theory  used  here  and  the  2-D  plane 
strain  model  are  very  close. 


3.1  ANALYTICAL  FORMULATION 


The  formulation  used  here  deals  with  the  following  three  aspects  : 

.  Analysis  of  the  plate.  The  aim  is  to  get  strains  and  stresses  in  any 
point  of  the  plate. 

.  Failure  criterion.  Once  the  stress  tensor  is  known,  a  failure 
criterion  is  needed  in  order  to  know  the  relationship  between  the 
stress  tensor  and  the  strength/stress  ratio. 


.  Optimization  procedure.  The  objective  is  to  find  the  minimum  weight 
plate,  using  strength  as  a  design  criterion. 


3.1.1  Equations  of  analysis 
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The  analysis  of  the  structure  has  been  carried  out  by  means  of  the  finite 
element  method.  By  means  of  this  technique,  strains  and  stresses  can  be 
obtained  in  any  point  of  the  plate. 

With  the  finite  element  method,  the  interior  of  the  plate  is  idealized  as  an 
assembly  of  discrete  elements  over  which  the  unknown  displacements  are 
represented  approximately  by  linear,  quadratic,  etc.,  variations.  Thus,  the 
governing  equation  of  internal  equilibrium  is  satisfied  approximately. 

Let  5ae  be  a  virtual  displacement  of  the  nodes.  This  displacement  generates 
the  following  displacements  inside  the  finite  element: 

8u  =  N  8ae  (3.1) 

and  the  following  strains  inside  the  finite  element: 

8e  =  B  8ae  (3.2) 

The  work  done  by  the  node  forces  is  equal  to  the  addition  of  the  products 
of  each  component  of  force  by  the  corresponding  displacement: 

8aeT  qe  (3.3) 

Analogously,  the  internal  work-normalized  to  volume-done  by  distributed 
stresses  and  forces  is 


8eTa-8uT 


(3.4) 


or 

SaT  (BTa  -  NTb)  (3.5) 

The  external  work  must  be  equal  to  the  internal  work  integrated  on  the 
elemental  volume.  Therefore, 
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(3.6) 


I  Bto  dv  -  j  N1 
/V  Jv 


8aeTqe  =  8aeT  (  Btg  dv  -  |  N‘b  dv  ) 

'v° 

This  expression  is  valid  for  any  virtual  displacement.  Thus, 


=  J  BTo  dv  -  J  N1 
Jv  Jv1 


b  dv 


(3.7) 


For  a  linear  relationship  between  stresses  and  strains,  the  following 
expression  can  be  obtained: 


qe=  keae  +  f 


where  ke  is  the  elemental  stiffness  matrix 


(3.8) 


e=  J  BtD 
7v 


B  dv 


(3.9) 


and 


=  -  I  NTb  dv  - 1  BTDeo  dv  +  I 
Jv  JV>  J\ 


B  Deodv+  |  B^odv 

V" 


(3.10) 


The  force-displacement  relation  for  the  whole  plate  is  given  by: 

[K]{a)=q  (3.11) 

where  [K]  is  the  stiffness  matrix  of  the  structure,  {a}  is  the  displacement 
vector,  and  {q}  is  the  force  vector. 

The  strain  energy  eij  of  the  jth  layer  of  the  ith  element  is  given  by: 

eij  =  {a}il  [kjy  {a)i  (3.12) 
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The  element  strains  {e}j  in  the  ith  element  are  expressed  by: 


{e}p  [c]  {a)i  (3.13) 

where  [c]  is  the  strain-displacement  matrix.  The  stresses  {a}  ij  in  the  jth 
layer  of  the  ith  element  are  given  by: 

{or} r  [Q]ij  {eh  (3.14) 

where  [Qhj  is  the  matrix  of  elastic  constants  of  the  jth  layer  of  the  ith 
element. 


3.1.2  Failure  criterion 


Once  we  know  how  to  calculate  strains  and  stresses  in  any  point  of  the 
plate,  we  need  to  apply  a  failure  criterion  in  order  to  assess  the  stress  level  in 
each  layer. 

For  optimization  processes,  the  concept  of  strength/stress  ratio  is 
extremely  useful.  Thus,  the  chosen  failure  criterion  must  have  a  unique 
strength/stress  ratio  for  each  combined  state  of  stress  and  the  corresponding 
state  of  strain.  For  this  reason,  a  quadratic  criterion  has  been  applied. 

According  to  Tsai  and  Wu  [72],  there  exists  a  failure  surface  in  the  stress- 
space  in  the  following  scalar  form: 

Fij  Oy  +  Fi  0i  =  1  i,j  =  1,2, 3, 4, 5, 6  (3.15) 

where  the  contracted  notation  is  used;  Fj  and  Fjj  are  strength  tensors  of  the 
second  and  fourth  rank,  respectively. 
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One  method  of  applying  the  strength  criterion  is  to  transform  the  stress 
components  into  the  material-symmetry  axes.  The  final  expression  will  be 
composed  by  the  following  terms: 

Fj  Oi  +  F2  ((?2  "*■  ^3  ^  +  Fn  G  1  2 

+  F22*  22  +°  3  2  +  2  G  42  )  +  F66  (g  52  +G  g'2) 

+  2  F12  (<3\  02'  +  Gi'  03' )+  2F23(g2'  G3  +  G  4  2)  =  1  (3.16) 

where 


Fi  =  1/X-l/X' 

(3.17) 

Fn  =  i/(X  x1) 

(3.18) 

F2  =  1/Y-l/Y’ 

(3.19) 

F22  =  1/(Y  Y') 

(3.20) 

F33  =  1/(Z  Z') 

(3.21) 

f66  -  1/(S  S’) 

(3.22) 

The  coupling  or  interaction  terms 

can  be  calculated  by  means  of  the 

following  formula: 

Fij  =  Fy  *  [  FU  Fjj  ]  1/2  i#j 

(3.23) 

and  Fy  *  =  -0.5 

(3.24) 

The  strength/stress  ratio  r  is  the  ratio  between  the  maximum,  ultimate  or 
allowable  strength,  and  the  applied  stress: 

{^Imax  =  r  applied  (3.25) 
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Since  each  combination  of  stress  components  in  Equation  3.26 


Fy  ai  Oj  +  Fi  ai  =1  (3.26) 

reaches  its  maximum  when  the  left-hand  side  reaches  unity,  we  can  substitute 
Equation  3.25  into  Equation  3.26  : 

[Fy  <Ji  Oj  ]  r2  +  [Fi  <jj]  r  - 1  =  0  (3.27) 

If  we  define 

a  =  Fy  Oj  Cj  (3.28) 

b  =  Fi  <Ji  (3.29) 

Then,  the  strength/stress  ratio  r  is  the  positive  square  root  in  the  quadratic 
formula  : 

r  =  -  (  b  /2a  )  +  [(  b  /2a  )2+  1/a  ]l/2  (3.30) 


3.1.3  Recurrence  relations 


The  recurrence  relations  proposed  here  for  resizing  the  elements  are  based 
on  the  optimality  criteria,  which  are  not  rigorous  in  the  mathematical  sense, 
but  which  are  found  to  give  near  optimum  weight  designs  for  large 
structures  in  a  very  efficient  way. 

The  optimality  criterion  for  the  generalized  stiffness  requirement  can  be 
stated  as,  "the  optimum  structure  is  the  one  in  which  the  ratio  of  the  average 
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strain  energy  density  to  the  mass  density  is  the  same  in  all  the  elements"  [73]. 
This  criterion  can  be  written  as: 

1  =  X  ey  /  pij  i=l  j _ ,  m 

j=l,_ ,  m  (3.31) 

where  X  is  the  Lagrangian  parameter,  ey  is  the  strain  energy  density  of  the 
jth  layer  of  the  ith  element,  and  py  is  the  mass  density.  The  strain  energy 
density  is  given  by: 


eij  —  eij  /Vy  (3.32) 

where  Vy  is  the  volume  of  the  element  defined  by: 

Vy  =  ty  li  (3.33) 

where  ty  is  the  thickness  of  the  jth  layer  of  the  ith  element  and  \  is  the  surface 
area  of  the  ith  element.  The  design  variable  ty  can  be  written  as: 

tjj  =  A  Oij  (3.34) 

where  ay  is  the  relative  thickness  (normalized  to  the  maximum  thickness)  of 
the  jth  layer  of  the  ith  element,  and  A  is  the  scaling  parameter.  Introducing 
the  scaling  parameter  in  Equation  3.11  gives: 

A  [K]  [a]  =  [q]  '  (3.35) 

or 

[K’l  {a’}  =  [q]  (3.36) 


23 


where 


{a}  =  {a'}  /  A  (3.37) 

and  [K’l  is  the  stiffness  matrix  for  the  whole  structure  obtained  by  using  the 
relative  design  vector  ay.  Introducing  the  scaling  parameter  intro  Equations 

3.12  to  3.14,  the  relations  between  the  actual  quantities  and  the  relative 


quantities  at  element  level  can  be  expressed  as: 

[Kljj  =  A  [K'lij  (3.38) 

{ a)  j  =  {a'}j/ A  (3.39) 

{eh=  {e'h/A  (3.40) 

{o)ij=  {o’Jij/ A  (3.41) 

where  the  prime  quantities  are  the  relative  values. 

Introducing  Equations  3.32,  3.38  and  3.39  in  Equation  3.31  gives: 

1  =  X  e  y  /  (  A2  py)  (3.42) 

where 

e*  ij  =  {r’Jtj  [K'lij  ( r* } i  /  (ay  li)  (3.43) 

Multiplying  both  sides  of  Equation  3.42  by  a2y  and  taking  the  square  root 
gives: 

ay  =  B  ay  (e’y  /  py)1/2  (3.44) 

where  B  is  a  constant.  Equation  3.44  can  be  rewritten  in  an  iterative  form  as: 
(a  ij  )v+l  =  B  (ct  jj  )v  (e  jj  /  p  yX'1^  (3.45) 
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where  v  and  v+1  refer  to  the  cycles  of  iteration. 


In  this  procedure  the  resizing  of  an  element  is  done  by  dividing  the 

design  variable  by  the  minimum  actual  strength/stress  ratio  for  that  element: 

* 

(aij)v+l  ~  (®ij)v  /  OOminv  (3.46) 

Therefore,  the  election  of  the  appropriate  failure  criterion  is  critical.  If 
the  strength/stress  ratio  defined  by  the  failure  criterion  is  unique  for  each 
combined  state  of  stress  and  the  corresponding  state  of  strain,  the  application 
of  Equation  3.46  is  efficient  and  simple.  Otherwise,  the  application  of  the 
procedure  becomes  more  complicated. 


3.2  ANALYSIS  OF  THE  VARIABLE  THICKNESS  PROBLEM 


The  aim  of  this  study  is  to  find  out  what  happens  inside  the  laminate  when 
the  laminate  thickness  varies  and  one  transverse  load  is  applied.  The  behavior 
of  a  two-dimensional  plate  with  variable  thickness  is  extremely  difficult  to 
understand  because  it  is  a  3-D  problem,  and  because  bending,  shear  and 
variable  thickness  effects  appear  simultaneously.  In  other  words,  the  whole 
stress  tensor  must  be  considered:  the  three  in-plane  stress  components  are 
present  owing  to  the  bending  effects,  and  so  are  the  two  interlaminar  shear 
stress  components  because  of  the  shear  effects.  Finally,  the  interlaminar 
normal  stress  component  must  also  be  considered  because  of  the  variable 
thickness  effect. 

In  order  to  simplify  the  problem  and  to  understand  the  variable  thickness 
phenomenon,  a  one-dimensional  plate  with  variable  thickness  is  studied  here 
by  means  of  a  2-D  plane  strain  finite  element  model. 

The  next  section  refers  to  the  attempts  for  obtaining  through-thickness 
stress  distributions  in  a  cross  section  where  the  laminate  thickness  varies. 
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Once  the  stress  level  is  known,  the  laminate  failure  will  be  predicted  by 
means  of  a  quadratic  failure  criterion.  The  aim  of  this  study  is  to  analyze 
theoretically  the  influence  of  the  different  parameters  on  the  mechanical 
behavior  of  variable-thickness  laminated  composite  plates. 

A  study  of  the  damage  in  variable  thickness  laminated  composites  caused 
by  a  transverse  load  provided  an  understanding  og  composite  failure 
mechanisms  (caused  by  transverse  loads)  and  identified  the  essential 
parameters  causing  the  damage  in  composite  plates. 

Finally,  an  experimental  study  will  assess  the  accuracy  of  the  theory,  as 
well  as  analyze  the  influence  of  the  different  parameters  on  the  mechanical 
behavior  of  variable-thickness  laminated  composite  plates. 

3.2.1  Model  assumptions  and  method  of  analysis 

The  model  used  here  assumes: 

.  This  study  is  limited  to  one-dimensional  laminated  composite  plates,  in 
order  to  avoid  three-dimensional  effects  and  to  analyze  properly  the 
consequences  of  the  variable  thickness  effect. 

.  The  fiber  orientation  is  longitudinal  (x-direction).  This  is  the  optimal 
direction  for  one-dimensional  laminated  composite  plates. 

.  The  material  used  is  graphite/epoxy  AS4/3501-6. 

The  method  of  analysis  is  explained  by  the  following  points: 

.  A  2-D  plane  strain  finite  element  model. 

.  The  mesh  is  composed  by  1891  nodes. 


.  Different  types  of  loading  and  boundary  conditions  have  been  studied 
to  analyze  the  variable  thickness  problem  from  a  theoretical  point  of 
view  as  well  as  to  obtain  through-thickness  stress  distributions. 

.  A  three-point  bending  model  has  been  used  for  analyzing  theoretically 
the  influence  of  the  different  parameters  on  the  mechanical  behavior  of 
variable-thickness  laminated  composite  plates. 

.  A  three-point  bending  testing  in  fatigue  and  static  conditions  have  been 
carried  out  for  identifying  the  failure  mechanisms  . 

.  A  three-point  bending  testing  in  static  conditions  has  been  used  to 
verify  the  theoretical  model,  and  to  analyze  experimentally  the 
influence  of  the  different  parameters  on  the  mechanical  behavior  of 
variable-thickness  laminated  composite  plates. 


3. 2. 2. Through-thickness  stresses  distributions 


This  section  explains  the  attempts  to  obtain  through-thickness  stresses 
distributions  in  a  cross  section,  where  the  laminate  thickness  varies  from  t  to 
t-ti. 

The  solution  of  the  problem  is  a  function  of  a  number  of  parameters  : 

.  0  :  angle  of  variation  of  thickness 
.  t/ti  :  thickness  ratio 
.  Stacking  sequence 
.  Bending  moment  distribution 
.  Shear  force  distribution 
.  Characteristic  length  /  laminate  thickness  ratio 
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.  Material. 


The  model  described  in  Section  3.2.1  has  been  applied  to  different  cases  in 
order  to  find  out  what  happens  inside  the  laminate  when  the  laminate 
thickness  varies  and  to  study  the  sensitivity  of  the  the  parameters  mentioned 
above.  Thus,  once  the  relationship  between  each  parameter  and  the 
mechanical  behavior  of  the  plate  is  known,  some  design  changes  can  be  made 
in  order  to  improve  the  mechanical  characteristics  of  the  plate. 

Figures  3,  4,  and  5  show  the  distributions  of  03  and  CJ5  through  the 
thickness  for  one-dimensional  two-sides  tapered  plates  and  linear  bending 
moment  distributions  with  9M/dx  >  0.  These  stress  components  present 
remarkable  differences  with  respect  to  the  constant-thickness  plate 
distributions.  Especially,  those  related  to  <73  and  05^  because  05  has  two 
maxima  near  the  top  and  bottom  surfaces  and  <73  reaches  much  higher  values 
than  in  the  case  of  constant  thickness. 

In  Figures  6,  7,  and  8,  the  distributions  of  <7^  and  <55  through  the 
thickness  for  one-dimensional  one-side  tapered  plates  and  linear  bending 
moment  distributions  with  9M/3x  >  0.  Obviously,  to  tape  the  plate  in  just  one 
side  leads  to  non-symmetrical  distributions  with  peaks  of  interlaminar  stresses 
near  the  surface  tapered. 

Figures  9,  10,  and  11  show  the  distributions  of  03^  and  05  through  the 
thickness  for  one-dimensional  two-sides  tapered  plates  and  linear  bending 
moment  distributions  with  3M/3x  <  0.  The  sign  of  9M/9x  has  a  strong 
influence  on  the  distributions  of  longitudinal  and  interlaminar  stresses. 

Finally,  Figures  12,  13,  and  14  represent  the  distributions  of  (7^  (73,  and  <75 
through  the  thickness  for  one-dimensional  one  side  tapered  plates  and  linear 
bending  moment  distributions  with  3M/3x  <  0. 

In  Figures  15,  16,  and  17,  <7it  <73,  and  05  distributions  are  represented  for 
t i/t=0.3,  M/Pt=5,  and  3M/3x  <  0,  respectively.  Two  angles  of  variation  of 
thickness  have  been  considered:  15  and  90  degrees.  To  reduce  the  angle  of 
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variation  of  thickness  from  90°  to  15°  leads  to  decrease  considerably  the 
values  of  the  peaks  in  the  distributions  of  <73  and  <75,  especially  these  last 

two. 


Figures  18,  19,  and  20,  show  the  distributions  of  <Ji,  <73,  and  <75  for 
ti/t=0.2,  M/Pt=5,  and  3M/3x  <  0,  respectively.  As  expected,  maxima  values 

decrease  considerably  with  respect  to  Figures  15,  16  and  17,  because  of  the 
smaller  value  of  ti/t. 

Distributions  of  <7i(  and  <75  for  ti/t=0.1,  M/Pt=5,  and  3M/3x  <  0  are 

presented  in  Figures  21,  22,  and  23,  respectively.  Numerical  values  keep 
decreasing  because  of  the  value  of  tj/t=0. 1 .  Also,  the  effect  of  the  angle  of 

variation  of  thickness  can  be  seen  by  comparing  the  distributions  for  0=90° 
and  0=15°. 

The  following  conclusions  can  be  drawn  after  analyzing  the  results 
described  above: 

.  The  examination  of  Figures  3  through  14  leads  to  the  conclusion  that  the 
variable-thickness  effect  presents  a  strong  influence  on  the  distributions  of 
longitudinal  and  interlaminar  stresses.  High  peaks  of  interlaminar  stresses 
appear  near  the  areas  of  change  of  thickness. 

.  Comparison  of  Figures  3,  4,  and  5  and  Figures  6,  7,  and  8  shows  that  the 
distributions  of  the  through-thickness  stresses  vary  in  function  of  the  sign 
of  3M/3x. 

.  Stress  distributions  are  strongly  dependent  on  t/tl  or  thickness  ratio,  as  the 
comparisons  of  Figures  15  through  23  reflect. 

.  0,  or  angle  of  variation  of  thickness  has  a  remarkable  influence  on 

interlaminar  stress  distributions  through  the  laminate  thickness.  This  fact 
is  shown  in  Figures  15  through  23. 
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Focusing  on  optimization  of  laminated  composite  plates,  the  last  four 
parameters  in  the  list  mentioned  at  the  beginning  of  section  3.2.2  will  be  input 
data,  and  the  optimization  procedure  will  determine  the  rest  of  parameters  of 
the  list  except  for  the  angle  of  variation  of  thickness.  Therefore,  the 
sensitivity  study  of  this  parameter  becomes  more  and  more  important  since 
the  angle  of  variation  of  thickness  is  critical  for  the  mechanical  behavior  of 
variable  thickness  plates. 


3.2.3  Failure  mechanisms 

A  variable-thickness  unidirectional  laminated  composite  exhibits  a  variety 
of  failure  mechanisms  when  subjected  to  transverse  loads.  The  knowledge  of 
these  modes  is  very  important  for  the  understanding  of  what  happens  inside 
the  laminate.  Once  we  know  which  the  critical  failure  mechanisms  are,  we 
will  be  able  to  identify  the  strain  components  associated  with  them. 

A  number  of  laminated  composite  plates  were  tested  in  fatigue  and  static 
conditions: 

.  AS4/3501-6  Graphite/Epoxy  was  used  for  this  work. 

.  The  types  of  testing  performed  were  three-point  static  and  dynamic 
bending. 

.  The  load  rate  for  the  static  testing  was  0.06  in/min  (0.001524 
m/min). 

.  For  the  fatigue  testing,  the  frequency  was  5  cycles/second,  and  the 
load  was  60  %  of  the  static  failure  load.  The  vertical  displacement 
was  controlled. 

Figure  24  shows  the  three-point-bending  testing. 
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After  studying  the  results,  the  following  failure  mechanisms  have  been 
reported: 

1 .  Delamination.  This  is  the  most  usual  failure  mode  when  the  angle  of 
variation  of  thickness  is  high.  It  is  due  to  a  sharp  thickness  variation  that 
generates  a  peak  of  interlaminar  shear  stress  near  the  areas  of  change  of 
thickness.  In  Figure  25,  a  delamination  failure  in  a  tapered  surface  is 
shown.  This  mechanism  can  also  be  detected  in  thick  plates,  though  the 
angle  of  variation  of  thickness  is  low.  A  representation  of  this  failure 
mechanism  can  be  seen  in  Figure  26.  Figure  27  shows  a  SEM  image  of  a 
delamination  fracture  surface.  This  damage  was  originated  by  a 
interlaminar  shear  stress  near  the  midplane  of  the  laminate.  A  detail  of  this 
fracture  surface  can  be  seen  in  Figure  28. 

2.  Bending  modes  (compression).  This  mechanism  occurs  when  the  plate  is 
very  thin,  the  variable  thickness  effect  is  negligible,  and  the  compression 
strength  is  lower  than  the  tension  strength  in  the  fiber  direction.  The 
following  in-plane  compression  failure  modes  can  be  distinguished: 

2.1.  Induced  transverse  tensile  failure.  Unidirectional  composites  can 
fracture  along  the  fibers  when  loaded  by  compression  by  a  transverse 
tensile  failure  mode,  because  of  the  weakness  of  the  matrix  and  the 
fiber-matrix  interface,  compared  with  the  strength  of  the  fibers. 

2.2.  Compressive  delamination  failure.  If  a  fiber  buckles,  the  fiber-matrix 
interface  may  fracture  in  shear  and  lead  to  ultimate  failure. 

2.3.  Euler  failure.  If  the  matrix  is  ductile  and  the  interface  is  strong,  the 
fiber  can  bend  without  matrix  failure  and  will  eventually  fracture  in 
bending. 

2.4.  Microbuckling.  A  more  likely  failure  mode  of  unidirectional  composite 
laminates  associated  with  fiber  microbuckling  and  fiber  kinking,  is 
shear  crippling.  Macroscopically,  shear  crippling  looks  like  a  shear 
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failure  on  a  plane  at  an  angle  to  the  direction  of  loading.  Microscopic 
inspection,  however,  indicates  that  shear  crippling  is  frequently  the 
result  of  kink-band  formation. 

2.5.  Strength  failure.  The  final  failure  mode  exhibited  in  unidirectional 

composites  is  associated  with  pure  compression  failure  of  the  fibers.  In 
this  case,  the  fracture  surface  is  likely  to  be  at  an  angle  to  the  loading 
direction,  usually  about  45°. 

3.  Bending  modes  (tension).  This  mechanism  occurs  when  the  plate  is  very 
thin,  the  variable  thickness  effect  is  negligible,  and  the  compression 
strength  is  higher  the  than  tension  strength  in  the  fiber  direction.  The 
following  in-plane  tension  failure  modes  can  be  distinguished: 

3.1.  Unidirectional  composite  subjected  to  longitudinal  tensile  load.  Brittle 
failure.  In  this  mode,  stress  concentrations  created  at  the  broken  fiber 
ends  will  lead  to  specimen  separation  at  a  given  cross  section. 

3.2.  Unidirectional  composite  subjected  to  longitudinal  tensile  load.  Brittle 
failure  with  fiber  pullout.  Variations  in  bond  strength  and  local  load 
transfer  mechanisms  from  matrix  to  fiber  can  lead  to  the  pull-out  of 
the  fibers  from  the  matrix  at  fracture. 

3.3.  Unidirectional  composite  subjected  to  longitudinal  tensile  load.  Brittle 
failure  with  debonding  and/or  matrix  failure.  Finally,  in  other  cases, 
cracks  at  different  cross  sections  of  the  laminate  may  join  together  at 
fracture  through  fiber-matrix  debonding  or  by  shear  failure  of  the 
matrix.  This  interfiber  matrix  shear  failure  and  fiber-matrix  debonding 
can  occur  either  independently  or  in  combination;  that  is,  portions  of 
the  failure  path  may  exhibit  debonding,  while  matrix  shear  failure  is 
evident  in  other  regions, 

4.  Mixed  modes.  A  mixed  failure  mechanism  has  been  detected  in  the  regions 
close  to  the  application  of  the  load.  There  exists  a  compressive 
delamination  because  of  bending  mode  (compression),  and  crushing 
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because  of  compressive  peeling  or  interlaminar  normal  stress.  This 
mechanism  was  observed  in  the  three-point-bending  fatigue  testing,  after 
50,000  cycles.  Figure  29  shows  a  representation  of  this  mixed  mode. 

The  SEM  photograph  of  Figure  30,  points  clearly  the  fracture  area.  In 
Figures  31  through  35,  some  SEM  images  show  the  different  mechanisms. 
At  the  top  (Figure  31),  a  crushing  failure  because  of  compressive 
interlaminar  normal  stress  is  shown.  In  Figures  32  through  35,  some 
surface  views  reflect  a  clear  compressive  delamination  failure  because  of 
the  bending  mode  (compression).  No  kinking  was  reported  in  the  SEM 
images.This  mixed  mode  does  not  appear  in  the  four-point-bending  testing 
and  it  can  also  be  avoided  by  placing  a  hard  rubber  pad  (60-80  durometer) 
between  the  load  point  and  the  specimen. 

After  analyzing  the  failure  mechanisms  described  above,  the  following 
conclusions  can  be  drawn: 

.  The  critical  failure  modes  are  related  to  two  strain  components:  Ei  or 
longitudinal  normal  strain  and  £5  or  interlaminar  shear  strain. 

.  The  bending  (compression  or  tension)  modes  are  likely  to  occur  when  the 
plate  is  thin  and  the  variable  thickness  effect  is  negligible. 

.  The  delamination  mechanisms  are  usual  in  cases  where  the  variable 
thickness  effect  is  critical  (sharp  change  of  thickness).  The  failure  appears 
in  the  tapered  surface.  This  mechanism  also  occurs  when  the  plate  is  thick, 
even  if  the  variable  thickness  is  negligible.  In  this  case,  the  failure  appears 
near  the  middle  plane  of  the  laminate. 

.  The  mode  related  to  interlaminar  normal  strains  (mixed  mode)  is  due  to 
stress  concentrations  in  the  areas  near  the  application  of  the  load.  This 
mechanism  can  be  easily  avoided  by  means  of  the  solutions  mentioned 
above. 
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3.2.4  Experimental  study  and  verification  of  the  theoretical  model 


An  experimental  study  assessed  the  accuracy  of  the  theoretical  model,  and 
the  influence  of  the  angle  of  variation  of  thickness  on  the  strength  of  a 
composite  plate: 

.  AS4/3501-6  Graphite/Epoxy  was  used  for  this  work. 

.  The  type  of  testing  performed  was  three-point  static  bending. 

.  The  load  rate  was  0.06  in/min  (0.001524  m/min). 

.  Three  types  of  specimens  were  tested.  Maximum,  minimum  thickness, 
and  span  were  the  same  for  the  three  specimens.  The  only  variable 
parameter  was  the  angle  of  variation  of  thickness: 

.  Specimen  OLCP_7001.  Angle  of  variation  of  thickness:  6° 

.  Specimen  OLCP_7002.  Angle  of  variation  of  thickness:  45° 

.  Specimen  OLCP_7003.  Angle  of  variation  of  thickness:  90° 

.  The  span  was  8  inches  (0.2032  meters). 

.  The  maximum  thickness  was  at  the  center  of  the  beam:  0.66  inches 
(0.0167  meters). 

.  The  minimum  thickness  was  at  the  end  of  the  beam:  0.28  inches 
(0.007  meters). 

.  Longitudinal  and  interlaminar  strains  were  measured  by  using  strain 
gauges.  The  collocation  of  the  gauges  is  reported  in  Table  3.1.  The 
election  of  the  points  was  carried  out  in  function  of  the  theoretical 
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analysis.  Strain  gauges  were  disposed  in  those  points  whose  strain  value 
was  critical,  according  to  the  results  from  the  finite  element  method 
applied  to  variable  thickness  laminated  composite  plates. 

Table  3.1  Coordinates  x  and  z  of  strain  gauges.  Values  are  expressed  in  inches 
(meters). 


POINT 

X 

Z 

1 

ma 

(0.0000) 

MB 

(0.00254) 

2 

0.0 

(0.0000) 

■a 

(0.00508) 

3 

1.0 

(0.0254) 

MB 

(0.00889) 

4 

2.0 

(0.0508) 

0.415 

(0.01054) 

5 

3.0 

(0.0762) 

0.600 

(0.01524) 

6 

4.0 

(0.1016) 

■SB 

(0.01295) 

7 

4.0 

(0.1016) 

MB 

(0.00889) 

8 

4.0 

(0.1016) 

0.000 

(0.00000) 

The  study  on  failure  mechanisms  shows  that  critical  strain  components  are 
the  longitudinal  strain  ei,  and  the  interlaminar  shear  strain  is  £5.  An  analysis 

of  both  components  for  specimens  described  OLCP_7001,  OLCP_7002,  and 
OLCP_7003  follows. 

The  scheme  of  the  specimen  OLCP_7001  is  depicted  in  Figure  36. 
Measurements  are  reported  in  Figure  37.  This  specimen  was  made  by 
tapering  the  plate  in  its  inner  part  (Figure  38)  and  the  angle  of  variation  of 
thickness  is  6°.  The  failure  load  was  4950  lb  (22000N).  The  failure  mode  is 
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shown  in  Figure  39.  The  upper  surface  does  not  present  any  failure,  in  spite 
of  being  tapered.  Instead  of  that,  a  crack  appears  near  the  middle  plane  at  the 
end  of  the  plate  where  interlaminar  shear  strains  are  maxima.  The  reason  for 
this  is  that  the  angle  of  variation  of  thickness  is  very  low.  A  scheme  of  the 
situation  of  the  gauges  is  depicted  in  Figure  40.  The  numerical  results  of  the 
specimen  OLCP_7001  are  presented  in  Tables  3.2,  3.3,  and  3.4.  Tables  3.2 
and  3.3  show  longitudinal  and  interlaminar  shear  strains  in  the  critical  points 
for  the  failure  load,  respectively  .  The  maximum  value  is  registered  in  point 
number  1  (£5=0.0131). 

Table  3.2  Longitudinal  strains  in  specimen  OLCP_7001.  Load:  4950  lb 
(22000  N). 


POINT 

€|  *  IE-3 

TESTING 

£|  *  IE-3 

THEORY 

8 

8.9 

8.6 

Table  3.3  Interlaminar  shear  strains  in  specimen  OLCP_7001.  Load:  4950  lb 
(22000  N). 


POINT 


£5  *  IE-3 
TESTING 

13.1  (FAILURE) 


9. 


€5  *  IE-3 
THEORY 

12.5 


8.1 


The  maximum  vertical  displacement  is  reported  in  Table  3.4.  In  both 
fields,  strains  and  displacements,  an  excellent  agreement  is  found  between  the 
data  and  the  prediction. 

Table  3.4  Maximum  vertical  displacement  in  specimen  OLCP_7001.  Load: 
4950  lb  (22000  N). 


POINT 

dz  TESTING 

d  z  THEORY 

8 

-0.173  (0.00439) 

-0.161  (0.00409) 

The  scheme  of  the  specimen  OLCP_7002  is  depicted  in  Figure  41. 
Measurements  are  reported  in  Figure  42.  The  angle  of  change  of  thickness  is 
45°.  Six  steps  were  designed  along  the  specimen.  The  failure  load  was  2960 
lb  (13150  N).  The  failure  mode  is  shown  in  Figure  43.  Owing  to  the  high 
angle  of  variation  of  thickness,  failure  appears  at  the  upper  surface  in  the 
thinnest  area  of  change  of  thickness.  Failure  mechanism  is  a  kind  of 
interlaminar  shear  mode.  As  theoretical  results  predict,  there  is  a  peak  of 
interlaminar  shear  strains  in  the  places  where  the  thickness  changes.  A  scheme 
of  the  situation  of  the  gauges  is  depicted  in  Figure  44.  The  numerical  results 
of  the  specimen  OLCP_7002  are  presented  in  Tables  3.5,  3.6,  and  3.7.  Tables 
3.5  and  3.6  show  longitudinal  and  interlaminar  shear  strains  in  the  critical 
points  for  the  failure  load,  respectively.  The  maximum  value  is  registered  in 
point  number  3  (65=0.0131). 


37 


Table  3.5  Longitudinal  strains  in  specimen  OLCP_7002.  Load:  2960  lb 
(13150  N). 


POINT 

£|  *  IE-3 

TESTING 

*  IE-3 

THEORY 

8 

4.4 

3.96 

Table  3.6  Interlaminar  shear  strains  in  specimen  OLCP_7002.  Load:  2960  lb 
(13150  N). 


POINT 

€5  *  IE-3 
TESTING 

€5  *  IE-3 
THEORY 

1 

6.0 

5.7 

3 

13.1  (FAILURE) 

12.5 

4 

11.0 

1 1.6 

5 

10.0 

9.8 

7 

2.8 

2.3 

The  maximum  vertical  displacement  is  reported  in  Table  3.7.  In  both 
fields,  strains  and  displacements,  an  excellent  agreement  is  found  between  the 
data  and  the  prediction. 
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Table  3.7  Maximum  vertical  displacement  in  specimen  OLCP_7002.  Load: 
2960  lb  (13150  N). 


POINT 

d2  TESTING 

dz  THEORY 

8 

-0.157  <0.004) 

-0.141  (0.0036) 

The  scheme  of  the  specimen  OLCP_7003  is  depicted  in  Figure  45. 
Measurements  are  reported  in  Figure  46.  The  angle  of  change  of  thickness  is 
90°.  Six  steps  were  designed  along  the  specimen.  The  failure  load  was  2485 
lb  (11046  N).  The  failure  mode  is  shown  in  Figure  47.  Owing  to  the  high 
angle  of  variation  of  thickness,  failure  appears  at  the  upper  surface  in  the 
thinnest  area  of  change  of  thickness  at  a  low  failure  load.  As  with  specimen 
OLCPJ7002,  failure  mechanism  is  a  kind  of  interlaminar  shear  mode.  As 
theoretical  results  predict,  there  is  a  remarkable  peak  of  interlaminar  shear 
strains  in  the  places  where  the  thickness  changes.  A  scheme  of  the  situation  of 
the  gauges  is  depicted  in  Figure  48.  The  numerical  results  of  the  specimen 
OLCP_7003  are  presented  in  Tables  3.8,  3.9,  and  3.10  .  Tables  3.8  and  3.9 
show  longitudinal  and  interlaminar  shear  strains  in  the  critical  points  for  the 
failure  load,  respectively.  The  maximum  value  is  registered  in  point  number 
3  (£5=0.0139). 


Table  3.8  Longitudinal  strains  in  specimen  OLCP_7003.  Load:  2485  lb 
(11046  N). 


POINT 


Z 


1  *  IE-3 
TESTING 


Z 


*  IE-3 
THEORY 


Table  3.9  Interlaminar  shear  strains  in  specimen  OLCP_7003.  Load:  2485  lb 
(11046  N) 


POINT 

£5  *  IE-3 
TESTING 

£5  *  IE-3 
THEORY 

1 

7.25 

7.5 

3 

13.9CFAILURE) 

12.5 

4 

12.18 

11.6 

5 

1 1.02 

9.8 

7 

3.77 

- —  -  — 

3.4 

The  maximum  vertical  displacement  is  reported  in  Table  3.7.  In  both 
fields,  strains  and  displacements,  an  excellent  agreement  is  found  between  the 
data  and  the  prediction. 

Table  3.10  Maximum  vertical  displacement  in  specimen  OLCP_7003.  Load: 
2485  lb  (11046  N). 


POINT 

dz  TESTING 

dz  THEORY 

8 

-0.171  (0.00432) 

-  0.157  (0.004) 

Tables  3.11,  3.12,  and  3.13  compare  the  three  specimens  tested.  These 
tables  show  the  values  of  longitudinal,  interlaminar  shear  strains  and  maxima 
vertical  displacements,  respectively,  for  a  load  of  1000  lb  (4445  N). 
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Table  3.11  Longitudinal  strains  in  specimens  OLCP_7001,  OLCPJ7002  and 
OLCPJ7003.  Load:  1000  lb  (4445  N). 


POINT 

*  IE-3 

0LCP—7001 

€,  *  IE-3 
0LCP—7002 

€,  *  IE-3 
0LCP-7003 

8 

1.78 

1.77 

1.76 

As  expected,  the  longitudinal  strains  at  the  inner  surface  in  the  middle  of 
the  span  are  very  similar  for  the  three  specimens  tested.  That  is  because  the 
three  specimens  present  the  same  thickness  at  the  middle  of  the  span  and, 
therefore,  the  longitudinal  strains  and  stresses  should  be  very  close  in  the 
three  cases. 

Table  3.12  Interlaminar  shear  strains  in  specimens  OLCP_7001,  OLCP_7002 
and  OLCP_7003.  Load:  1000  lb  (4445  N). 


POINT 

£5  *  IE-3 
0LCP_700 1 

€5  *  IE-3 
OLCP_7002 

€5  *  IE-3 
0LCP—7003 

1 

2.62 

2.5 

2.4 

3 

1.26 

4.8 

5.24 

4 

1.14 

4.2 

4.4 

5 

0.84 

3.8 

4.0 

7 

2.54 

1.3 

1.12 

Table  3.12  shows  very  nice  results.  On  the  one  hand,  the  maxima  values  are 
registered  in  specimens  OLCP_7002  and  OLCP_7003,  whose  angles  of 
variation  of  thickness  present  high  values  (45  and  90  degrees,  respectively). 
On  the  other  hand,  in  these  two  cases,  the  values  given  by  gauges  number  3, 
4,  and  5  are  critical,  which  means  that  the  upper  surface  is  the  critical  one 
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from  the  failure  point  of  view.  Instead  of  that,  specimen  OLCP_7001  does 
not  present  its  maximum  value  in  any  of  these  points  and,  therefore,  the 
tapered  surface  is  not  critical  at  all. 

Table  3.13  Maximum  vertical  displacement  in  specimens  OLCP_7001, 
OLCPJ7002  and  OLCP_7003.  Load:  1000  lb  (4445  N). 


POINT 

dz 

0LCP—7001 

d2 

0LCP-7002 

dz 

0LCP—7003 

8 

-0.0346  (0.00087) 

-  0.059  (0.00149) 

-0.0628(0.00159) 

The  results  reported  in  Table  3.13  point  out  a  higher  stiffness  of  specimen 
OLCP_7001  with  respect  to  the  other  two.  The  lower  the  angle  of  variation 
of  thickness,  the  better  the  fiber  works.  Hence,  according  to  these  results, 
there  is  no  linearity,  and  the  angle  of  variation  of  thickness  must  be  very  low 
in  order  to  get  a  high  quality  design  from  the  point  of  view  of  stiffness  and 
strength. 

Figures  50  and  51  compare  the  distributions  of  interlaminar  strains 
through  the  laminate  thickness.  Both  graphics  refer  to  the  three-point-bending 
testing  and  the  specimen  OLCP_7001  described  above.  Figure  50  shows  the 
theoretical  and  experimental  distributions  of  interlaminar  shear  strain  £5  in 

section  AA.  As  shown  in  Figure  49,  section  AA  corresponds  to  the  end  of  the 
beam.  This  section  is  critical  because  the  laminate  thickness  is  minimum. 
Figure  51  represents  the  same  distributions  of  £5  in  section  BB,  at  the  middle 

of  the  span  (Figure  50).  In  both  cases,  theoretical  and  experimental 
distributions  are  very  close. 

Comparison  of  these  results  shows: 

.  For  low  values  of  the  angle  of  variation  of  thickness,  the  strength  is 
outstanding  and  the  variable  thickness  effect  does  not  alter  the 
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mechanical  behavior  of  the  plate.  Failure  occurs  at  the  middle  plane  of 
the  plate  and  corresponds  to  a  delamination  mechanism  (owing  to 
interlaminar  shear  stresses). 


.  For  high  values  of  the  angle  of  variation  of  thickness,  the  strength  is 
low  and  the  variable  thickness  effect  does  alter  the  mechanical  behavior 
of  the  plate.  Failure  occurs  at  the  area  where  the  thickness  varies,  and 
corresponds  to  a  delamination  mechanism  (owing  to  high  interlaminar 
shear  stresses  provoked  by  the  variable  thickness  effect). 

.  The  numerical  results  from  theory  and  testing  are  very  close  both  in 
the  fields  of  deflections  and  longitudinal  and  interlaminar  strains. 

.  Theoretical  and  experimental  distributions  of  through-thickness  strains 
are  also  very  close. 


3.3  DESCRIPTION  OF  RESULTS 


This  section  describes  a  number  of  numerical  results  for  one-  and  two- 
dimensional  laminated  composite  plates.  First,  the  method  of  analysis.  The 
conclusions  drawn  in  the  last  section  are  the  key  for  the  formulation  of  the 
present  problem.  The  results  obtained  from  the  2-D  plane  strain  model  are 
very  close  to  the  data  given  by  the  experimental  study.  This  model  can  be 
directly  applied  for  analyzing  one-dimensional  laminated  composite  plates. 
And  that  model  can  also  be  the  base  for  a  more  general  one,  applicable  to  the 
analysis  of  two-dimensional  laminated  composite  plates. 

Once  the  method  of  analysis  is  known,  the  different  input  data  required 
for  the  calculation  must  be  obtained.  In  a  structural  optimization  problem, 
there  are  two  types  of  input  data:  geometric  and  those  related  to  the  material 
itself.  There  is  no  problem  in  getting  to  know  the  geometric  data  of  the  plate, 
but  the  elastic  constants  and  strengths  of  the  material  require  a  special 
treatment.  In-plane  constants  are  given  in  the  literature,  though  some  doubts 
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exist  about  compression  strength.  However,  through-thickness  properties  are 
hard  to  get,  and  even  some  constants  like  interlaminar  normal  and 
interlaminar  shear  moduli  are  a  function  of  the  stacking  sequence.  Therefore, 
the  elastic  constants  and  strengths  will  be  obtained  for  the  material  used  in 
this  work. 

Finally,  one-  and  two-dimensional  laminated  composite  plates  will  be 
studied,  and  a  number  of  figures  will  show  the  optima  configurations  for 
different  kinds  of  plates,  types  of  loading,  and  boundary  conditions. 


3.3.1  Model  assumptions  and  method  of  analysis 


In  the  last  section,  a  2-D  plane  strain  model  has  been  verified  by  means  of 
an  experimental  study.  Hence,  there  are  two  options  to  optimize  two- 
dimensional  plates: 

.  Generalizing  to  a  3-D  model,  by  using  a  three-dimensional  finite  element 
theory  and 

.  Using  a  2-D  model,  by  applying  a  shear  deformation  plate  theory. 

The  first  option  is  expensive,  especially  from  the  point  of  view  of  the 
optimization  and,  therefore,  the  subsequent  application  of  an  iterative 
procedure. 

The  second  option  is  efficient,  because  the  analysis  is  carried  out  very  fast 
and  the  required  number  of  nodes  of  the  mesh  is  not  large.  But  the  following 
question  remains  to  be  answered:  is  a  2-D  shear  deformation  plate  theory 
accurate  enough  to  study  two-dimensional  plates  with  variable  thickness  ? 
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This  question  has  two  possible  answers: 

.  If  the  angle  of  variation  of  thickness  is  90  degrees,  there  are 
discontinuities  in  the  exterior  surface(s)  of  the  laminated  composite  plates, 
there  are  free  edges  effects,  and  a  3-D  finite  element  theory  should  be 
used. 

.  If  the  angle  of  variation  of  thickness  is  low,  and  the  plate  exterior 
surfaces  are  continuous  as  shown  in  Figure  36,  a  two-dimensional  plate 
with  variable  thickness  can  be  studied  very  accurately  by  means  of  a  2-D 
shear  deformation  theory.  This  fact  is  based  on  a  number  of  verifications 
made  between  2-D  plane  strain  models  and  1  -D  models  using  a  shear 
deformation  plate  theory.  In  all  cases,  deflections  were  very  close  and, 
according  to  the  conclusions  of  failure  mechanisms,  the  critical  stress 
components  <3\  and  05  presented  very  small  differences.  One  verification 

is  presented  in  Figures  52,  53,  and  54.  The  geometry  is  shown  in  Figure 
52.  In  Figures  53  and  54,  distributions  of  Gj  and  05  from  both  models 

are  depicted,  respectively  . 

Both  theoretical  and  experimental  studies  about  variable  thickness  ( 
Section  3.2)  showed  that  the  angle  of  variation  of  thickness  has  a  strong 
influence  on  the  failure  load,  and  that,  the  lower  the  angle,  the  higher  the 
failure  load.  Hence,  it  is  logical  to  design  tapered  plates  with  low  angles  of 
variation  of  thickness  and  continuous  exterior  surfaces.  In  this  case,  a  2-D 
shear  deformation  plate  theory  is  applicable. 

The  finite  element  used  [67]  here  is  based  on  the  higher  order  shear 
theory  [62]  and  on  the  penalty  function  theory.  This  scheme  makes  it  possible 
to  analyze  thin  and  thick  plates,  because  of  its  general  formulation. 

The  results  given  by  the  finite  element  method  are  optimized  by  means  of 
the  iterative  procedure  described  in  Section  3.1.3.  Laminate  thickness  and 
fiber  orientation  are  the  design  variables.  Laminate  thickness  is  modified  in 
each  step  by  means  of  an  iterative  procedure  based  on  optimality  criteria  to 
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design  a  minimum  weight  structure.  In  order  to  assess  the  stress  level,  a 
quadratic  failure  criterion  is  applied  in  each  element.  If  the  stresses  in  all  the 
elements  satisfy  this  criterion,  the  process  is  over.  If  they  do  not,  another 
iteration  starts.  Usually,  around  ten  iterations  are  needed  to  get  convergence. 

The  optimization  of  fiber  orientations  is  complex  because  of: 

.  Practical  considerations.  Only  four  angles  have  been  used  : 

Oo,  450,  -45°*  and  90°  . 

.  Practical  considerations.  The  optimum  laminate  should  be  defined  in 
function  of  repetitive  sublaminates. 

.  Stacking  sequence.  A  strong  influence  on  the  strain  and  stress  level  is 
due  to  bending  and  interlaminar  effects. 

No  acceptable  result  was  obtained,  by  using  the  optimization  procedures 
available.  Hence,  2-D  shear  deformation  model  being  so  fast  and  playing 
with  discrete  variables,  all  the  possible  combinations  were  tried  as  possible 
optima  sublaminates. 

Finally,  the  following  assumptions  have  been  made  in  the  development  of 
the  analysis  : 

.  The  analysis  is  static. 

.  The  theory  applied  is  valid  for  small  deformations. 

.  The  material  used  is  supposed  to  be  elastic. 


3.3.2  Determination  of  elastic  constants  and  strengths 


46 


The  analysis  and  optimization  of  a  composite  structure  require  the 
determination  of  elastic  constants  and  strengths. 

T300/N5208  is  a  unidirectional  graphite/epoxy  laminate.  Thus,  each  layer 
can  be  considered  as  a  transversely  isotropic  material.  The  stiffness  matrix 
for  such  material  is  represented  in  Table  3.7.  There  are  five  independent 
constants  and  12  nonzero  components. 


Table  3.7  Stiffness  matrix  for  a  transversely  isotropic  material. 


£l 

£2 

£3 

£4 

£5 

£6 

C11 

C12 

C12 

0 

0 

0 

o2 

C21 

C22 

C23 

0 

0 

0 

c3 

C21 

C32 

C22 

0 

0 

0 

04 

0 

0 

0  (C22  -  C23)/2  0 

0 

0 

0 

0 

0 

c66 

0 

°6 

0 

0 

0 

0 

0 

C66 

The  five  independent  constants  present  the  following  values  : 


C11 

=  181.8 

GPa 

C12 

=  2.90 

GPa 

C22 

=  10.35 

GPa 

C23 

=  7.05 

GPa 

C66 

=  7.17 

GPa 

For  a  laminate,  the  components  that  relate  03  /  83  and  o 5  /es  are  a 
function  of  the  stacking  sequence.  According  to  Roy  and  Tsai  [74],  Figure  55 
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gives  some  values  of  the  interlaminar  normal  modulus  in  function  of  m  and  n 
in  a  laminate  T300/N5208  [90m  /0n]  s. 

Figure  56  [74]  shows  the  variation  of  the  interlaminar  shear  modulus  in 
function  of  <)>  in  a  laminate  T300/N5208  [+<j>  /-<}>]  s. 

The  strengths  for  this  material  are  given  by  the  following  values: 


X  = 

1500 

MPa 

X'  = 

1500 

MPa 

Y  = 

40 

MPa 

Y'  = 

246 

MPa 

Z  = 

40 

MPa 

Z'  = 

246 

MPa 

S  = 

68 

MPa 

S’  = 

68 

MPa 

3.3.3  One-dimensional  laminated  composite  plates 
The  following  points  analyze  the  problem: 

.  The  objective  is  to  get  the  minimum  weight  structure  by  using  strength  as 
a  design  criterion. 

.  The  laminate  used  is  [0]n  because  it  is  the  optimum  for  one-dimensional 
laminated  composite  plates  subjected  to  transverse  loads. 

.  Two-sided  tapered  laminates  have  been  considered. 

.  Three  kinds  of  plates  have  been  analyzed: 

-  thin  plates  (1/t  »  10), 

-  thick  plates  (1/t  «  10), 

-  sandwich  panels. 
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.  Boundary  conditions  are  restricted  to  : 

-  simply  supported  plates, 

-  clamped  plates, 

-  cantilever  plates. 

.  Two  types  of  loading  have  been  applied: 

-  uniform  load,  and 

-  point  load,  at  the  center  of  the  plate  for  simply  supported  and 
clamped  plates,  and  at  the  free  end  of  the  plate  for  cantilever  plates. 

.  The  mesh  was  composed  by  400  nodes. 

.  The  material  used  is  T300/N5208. 

Figure  57  shows  the  definition  of  coordinates  axis  and  the  key  to  colors  for 
one-dimensional  plates. 

Thin  plates  subjected  to  a  uniform  load 

Optimal  thin  plates  subjected  to  a  uniform  load  are  represented  in  Figure 
58.  Whatever  the  boundary  condition,  the  optimum  laminate  thickness  is  : 

t  =  2  k  a  (P/X)l/2  (3.47) 


where  a  is  the  plate  span,  P  is  the  load  applied,  X  is  the  uniaxial  tensile 
strength  of  a  ply  along  the  x-axis,  and  k  is  a  coefficient  that  varies  in  function 
of  x.  The  value  of  the  non-dimensional  coefficient  k  in  function  of  x  is  given 
for  the  following  boundary  conditions:  simply  supported,  clamped,  and 
cantilever  one-dimensional  plate.  Also,  the  weight  savings  (WS)  are  reported 
for  these  three  cases.  The  45.7  %  WS  for  the  cantilever  plate  is  remarkable. 
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Thick  plates  subjected  to  a  uniform  load 


Figure  59  shows  the  optimal  thick  plates  subjected  to  a  uniform  load.  In 
this  case,  the  optimum  laminate  thickness  is  : 

t  =  2  k  a  P/S  (3.48) 

where  S  is  the  shear  strength  in  the  xy-  or  12-plane  of  a  ply.  The  value  of  the 
non-dimensional  coefficient  k  in  function  of  x  and  the  optimal  configurations 
are  given  for  the  boundary  conditions  mentioned  above.  The  weight  saving  is 
the  same  for  the  three  boundary  conditions:  45.7  %. 

Sandwich  panels  subjected  to  a  uniform  load 

Optimal  one-dimensional  uniform  loaded  sandwich  panels  are  represented 
in  Figure  60.  In  this  case,  only  failure  mode  dealing  with  the  faces  have  been 
considered.  The  optimum  laminate  thickness  can  be  calculated  by  means  of 
the  following  expression: 

t  =  2  k  (a2  P/c  X)  (3.49) 

where  c  is  the  core  thickness.  Optimal  sandwiches  and  the  value  of  k  as  a 
function  of  x  are  presented  for  different  boundary  conditions.  The  weight 
saving  is  31.2  %  for  simply  supported  plates,  58.5  %  for  clamped  plates,  and 
56.9  %  for  cantilever  plates. 

Thin  plates  subjected  to  a  point  load 

Optimal  thin  plates  subjected  to  a  point  load  are  represented  in  Figure  61. 
Whatever  the  boundary  condition,  the  optimum  laminate  thickness  is  : 

t  =  2  k  (Pa/X)l/2  (3.50) 

The  value  of  the  non-dimensional  coefficient  k  in  function  of  x  is  given 
for  the  following  boundary  conditions:  simply  supported,  clamped,  and 
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cantilever  one-dimensional  plate.  Also,  the  weight  savings  (WS)  are  reported 
for  these  three  cases.  The  44.0  %  WS  for  clamped  plates  is  remarkable. 

Thick  plates  subjected  to  a  point  load 

Figure  62  shows  the  optimal  thick  plates  subjected  to  a  point  load.  In  this 
case,  the  optimum  laminate  thickness  is  : 

t  =  2  k  P/S  (3.51) 

The  value  of  the  non-dimensional  coefficient  k  in  function  of  x  and  the 
optimal  configurations  are  given  for  the  boundary  conditions  mentioned 
above.  There  is  not  weight  saving  in  this  case  for  any  boundary  condition. 

Sandwich  panels  subjected  to  a  point  load 

Optimal  one-dimensional  sandwiches  subjected  to  a  point  load,  are 
represented  in  Figure  63.  In  this  case,  only  failure  mode  dealing  with  the 
faces  have  been  considered.  The  optimum  laminate  thickness  can  be  calculated 
by  means  of  the  following  expression: 

t  =  2  k  (a  P/c  X)  (3.52) 

Optimal  sandwiches  and  the  value  of  k  as  a  function  of  x  are  presented  for 
different  boundary  conditions.  The  weight  saving  is  45.7  %  for  simply 
suppor  1  plates,  41.4  %  for  clamped  plates,  and  45.7  %  for  cantilever  plates. 

3.3.4  Two-dimensional  laminated  composite  plates 
The  following  points  analyze  the  problem: 

.  The  objective  is  to  get  the  minimum  weight  structure  by  using  strength  as 
a  design  criterion. 
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.  Only  symmetric  and  balanced  laminates  have  been  considered. 

.  Two-sided  tapered  laminates  have  been  considered. 

.  Three  kind  of  plates  have  been  analyzed: 

-  thin  plates  (1/t »  10), 

-  thick  plates  (1/t «  10), 

-  sandwich  panels. 

.  Boundary  conditions  are  restricted  to  : 

-  simply  supported  plates  along  the  four  edges  of  the  plate, 

-  clamped  plates  along  the  four  edges  of  the  plate. 

.  Two  types  of  loading  have  been  applied: 

-  uniform  load,  and 

-  point  load  at  the  center  of  the  plate. 

.  Owing  to  practical  considerations,  only  four  angles  have  been  used  : 
00,450^450,  and  900  . 

.  Plate  aspect  ratios  between  1  and  10  have  been  considered. 

.  The  mesh  was  composed  by  400  nodes. 

.  The  material  used  is  T300/N5208. 

Figure  64  shows  the  definition  of  coordinates  axis  and  key  to  colors  for 
one-dimensional  plates. 

Thin,  simply  supported  plates  subjected  to  a  uniform  load 


Weight  saving  and  normalized  deflection  for  thin,  simply  supported, 
uniform  loaded  plates  are  shown  in  Figure  65.  The  weight  saving  is  a 
function  of  the  aspect  ratio,  the  maximum  value  being  27%  for  a  square  plate. 
The  optimum  sublaminate  also  varies  in  the  function  of  the  aspect  ratio: 

.  [45/-45]  for  b/a  between  1  and  1.75, 

.  [0/45/-45]  for  b/a  between  1.75  and  2.25, 

.  [05/45/-45]  for  b/a  between  2.25  and  3.5, 

.  [0]  for  b/a  higher  than  3.5. 

The  quasi-isotropic  sublaminate  [0/45/-45/90]  presents  low-weight  saving 
with  respect  to  the  rest  of  sublaminates,  maximum  value  being  18  %  for  a 
square  plate. 

The  optimum  laminate  thickness  can  be  calculated  by  means  of  the 
following  expression: 

t  =  2  k  a(P/X)l/2  (3.53) 

where  a  is  the  length  of  the  plate,  P  is  the  load  applied,  X  is  the  uniaxial 
tensile  strength  of  a  ply  along  the  x-axis,  and  k  is  a  non-dimensional 
coefficient  that  is  a  function  of  x  and  y.  Figures  66  through  69  show  the 
values  of  k  along  the  x-axis  and  the  diagonal. 

Thin,  clamped  plates  subjected  to  a  uniform  load 

Weight  saving  and  normalized  deflection  for  thin,  clamped,  uniform 
loaded  plates  are  shown  in  Figure  70.  The  weight  saving  is  a  function  of  the 
aspect  ratio,  the  maximum  value  being  48  %  for  b/a=2.  The  optimum 
sublaminate  also  varies  in  the  function  of  the  aspect  ratio: 

.  [0/90]  for  b/a  between  1  and  1.5, 

.  [O5/9O2]  for  b/a  between  1.5  and  4.4, 
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.  [0]  for  b/a  higher  than  4.4. 

The  quasi-isotropic  sublaminate  [0/45/-45/90]  presents  low-weight  saving 
with  respect  to  the  rest  of  sublaminates,  maximum  value  being  40  %  for  a 
square  plate. 

The  optimum  laminate  thickness  can  be  calculated  by  means  of  the 
following  expression: 

t  =  2  k  a  (P/X)l/2  (3.54) 

Figures  71,  72,  and  73  show  the  values  of  the  non-dimensional  parameter  k 
along  the  x-axis  and  the  diagonal. 


Thick,  simply  supported  plates  subjected  to  a  uniform  load 

Weight  saving  and  normalized  deflection  for  thick,  simply  supported, 
uniform  loaded  plates  are  shown  in  Figure  74.  The  weight  saving  is  a 
function  of  the  aspect  ratio,  the  maximum  value  being  46  %  for  b/a=10.  The 
optimum  sublaminate  also  varies  in  the  function  of  the  aspect  ratio: 

.  [45/-45]  for  b/a  between  1  and  3.5, 

.  [0]  for  b/a  higher  than  3.5. 

The  quasi-isotropic  sublaminate  [0/45/-45/90]  presents  low-weight  saving 
with  respect  to  the  rest  of  sublaminates,  the  maximum  value  being  20  %  for 
b/a=10. 

The  optimum  laminate  thickness  can  be  calculated  by  means  of  the 
following  expression: 

t  =  2  k  a  P/S  (3.55) 
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Figures  75  and  76  show  the  values  of  the  non-dimensional  parameter  k 
along  the  x-axis  and  the  diagonal. 

Thick ,  clamped  plates  subjected  to  a  uniform  load 

Weight  saving  and  normalized  deflection  for  thick,  clamped,  uniform 
loaded  plates  are  shown  in  Figure  77.  The  weight  saving  is  a  function  of  the 
aspect  ratio,  the  maximum  value  being  46  %  for  b/a=10.  The  optimum 
sublaminate  also  varies  in  the  function  of  the  aspect  ratio: 

.  [45/-45]  for  b/a  between  1  and  3.5, 

.  [0]  for  b/a  higher  than  3.5. 

The  quasi-isotropic  sublaminate  [0/45/-45/90]  presents  low-weight  saving 
with  respect  to  the  rest  of  sublaminates,  the  maximum  value  being  20  %  for 
b/a=10. 

The  optimum  laminate  thickness  can  be  calculated  by  means  of  the 
following  expression: 

t  =  2k  a  P/S  (3.56) 


Figures  78  and  79  show  the  values  of  the  non-dimensional  parameter  k 
along  the  x-axis  and  the  diagonal. 

Simply  supported  sandwich  panels  subjected  to  a  uniform  load 

Weight  saving  and  normalized  deflection  for  simply  supported,  uniform 
loaded  sandwiches  are  shown  in  Figure  80.  The  weight  saving  is  a  function 
of  the  aspect  ratio,  the  maximum  value  being  35  %  for  b/a=10.  The  optimum 
sublaminate  also  varies  in  the  function  of  the  aspect  ratio: 

.  [45/-45]  for  b/a  between  1  and  2, 
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.  [0]  for  b/a  higher  than  2. 


The  quasi-isotropic  sublaminate  [0/45/-45/90]  presents  low  weight  saving 
with  respect  to  the  rest  of  sublaminates,  the  maximum  value  being  20  %  for 
b/a=10. 

The  optimum  laminate  thickness  can  be  calculated  by  means  of  the 
following  expression: 

t  =  2  k  (a2  P/c  X)  (3.57) 

Figures  81  and  82  show  the  values  of  the  non-dimensional  parameter  k 
along  the  x-axis  and  the  diagonal. 


Clamped  sandwich  panels  subjected  to  a  uniform  load 

Weight  saving  and  normalized  deflection  for  clamped,  uniform  loaded 
sandwiches  are  shown  in  Figure  83.  The  weight  saving  is  a  function  of  the 
aspect  ratio,  the  maximum  value  being  60  %  for  b/a=2.  The  optimum 
sublaminate  also  varies  in  the  function  of  the  aspect  ratio: 

.  [0/90]  for  b/a  between  1  and  1.5, 

.  [O5/9O2]  for  b/a  between  1.5  and  2.5, 

.  [0]  for  b/a  higher  than  2.5. 

The  quasi-isotropic  sublaminate  [0/45/-45/90]  presents  low-weight  saving 
with  respect  to  the  rest  of  sublaminates,  maximum  value  being  50  %  for 
b/a=10. 

The  optimum  laminate  thickness  can  be  calculated  by  means  of  the 
following  expression: 

t  =  2  k  (a2  P/c  X)  (3.58) 
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Figures  84,  85,  and  86  show  the  values  of  the  non-dimensional  parameter  k 
along  the  x-axis  and  the  diagonal. 

Thin,  simply  supported  plates  subjected  to  a  point  load 

Weight  saving  for  thin,  simply  supported  plates  subjected  to  a  point  load, 
is  shown  in  Figure  87.  The  weight  saving  is  a  function  of  the  aspect  ratio,  the 
maximum  value  being  76  %  for  b/a=2.  The  optimum  sublaminate  is  [0]  for 
all  plate  aspect  ratios  analyzed. 

The  quasi-isotropic  sublaminate  [0/45/-45/90]  presents  low-weight  saving 
with  respect  to  the  optimum  sublaminate,  the  maximum  value  being  63  %  for 
b/a=2. 

Since  no  general  formula  for  expressing  the  optimum  laminate  thickness 
has  been  found,  the  following  particular  case  has  been  calculated: 

.  a=7.87  in  (0.2  m), 

.  P=1.125  E4  lb  (5  E4  N). 

Figures  88,  89,  and  90  show  the  values  of  the  thickness  in  meters  along  the 
x-axis  and  the  diagonal. 

Thin,  clamped  plates  subjected  to  a  point  load 

Weight  saving  for  thin,  clamped  plates  subjected  to  a  point  load,  is  shown 
in  Figure  91.  The  weight  saving  is  a  function  of  the  aspect  ratio,  the 
maximum  value  being  83  %  for  b/a=10.  The  optimum  sublaminate  also  varies 
in  the  function  of  the  aspect  ratio: 

.  [45/-45]  for  b/a  between  1  and  1.5, 

.  [0]  for  b/a  higher  than  1 .5. 
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The  quasi-isotropic  sublaminate  [0/45/-45/90]  presents  low-weight  saving 
with  respect  to  the  optimum  sublaminate,  the  maximum  value  being  80  %  for 
b/a=10. 

Since  no  general  formula  for  expressing  the  optimum  laminate  thickness 
has  been  found,  the  following  particular  case  has  been  calculated: 

.  a=7.87  in  (0.2  m), 

.  P=1.125  E4  lb  (5  E4  N). 

Figures  92,  93,  and  94  show  the  values  of  the  thickness  in  meters  along  the 
x-axis  and  the  diagonal. 


Thick,  simply  supported  plates  subjected  to  a  point  load 

Weight  saving  for  thick,  simply  supported  plates  subjected  to  a  point  load, 
is  shown  in  Figure  95.  The  weight  saving  is  a  function  of  the  aspect  ratio,  the 
maximum  value  being  82  %  for  b/a=2  The  optimum  sublaminate  also  varies 
in  the  function  of  the  aspect  ratio: 

.  [45/-45]  for  b/a  between  1  and  2.25, 

.  [0]  for  b/a  higher  than  2.25. 

The  quasi-isotropic  sublaminate  [0/45/-45/90]  presents  low-weight  saving 
with  respect  to  the  optimum  sublaminate,  the  maximum  value  being  78  %  for 
b/a=2. 

Since  no  general  formula  for  expressing  the  optimum  laminate  thickness 
has  been  found,  the  following  particular  case  has  been  calculated: 

-  a=7.87  in  (0.2  m), 
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.  P=1.125  E4  lb  (5  E4  N). 


Figures  96  and  97  show  the  values  of  the  thickness  in  meters  along  the  x- 
axis  and  the  diagonal. 

Thick,  clamped  plates  subjected  to  a  point  load 

Weight  saving  for  thick,  clamped  plates  subjected  to  a  point  load,  is  shown 
in  Figure  98.  The  weight  saving  is  a  function  of  the  aspect  ratio,  the 
maximum  value  being  82  %  for  b/a=10.  The  optimum  sublaminate  also  varies 
in  the  function  of  the  aspect  ratio: 

.  [45/-45]  for  b/a  between  1  and  2.5, 

.  [0]  for  b/a  higher  than  2.5. 

The  quasi-isotropic  sublaminate  [0/45/-45/90]  presents  low-weight  saving 
with  respect  to  the  optimum  sublaminate,  the  maximum  value  being  78  %  for 
b/a=10. 

Since  no  general  formula  for  expressing  the  optimum  laminate  thickness 
has  been  found,  the  following  particular  case  has  been  calculated: 

.  a=7.87  in  (0.2  m), 

.  P=1.125  E4  lb  (5  E4  N). 

Figures  99  through  101  show  the  values  of  the  thickness  in  meters  along  the 
x-axis  and  the  diagonal. 

Simply  supported  sandwich  panels  subjected  to  a  point  load 

Weight  saving  for  simply  supported  sandwiches  subjected  to  a  point  load, 
is  shown  in  Figure  102.  The  weight  saving  is  a  function  of  the  aspect  ratio, 
the  maximum  value  being  94  %  for  b/a=2.  The  optimum  sublaminate  is  [0] 
for  all  the  plate  aspect  ratios  analyzed. 


59 


The  quasi-isotropic  sublaminate  [0/45/-45/90]  presents  low-weight  saving 
with  respect  to  the  optimum  sublaminate,  maximum  value  being  80  %  for 
b/a=2. 

Since  no  general  formula  for  expressing  the  optimum  laminate  thickness 
has  been  found,  the  following  particular  case  has  been  calculated: 

.  a=7.87  in  (0.2  m), 

.  P=1.125  E4  lb  (5  E4  N). 

In  Figures  103  through  105  show  the  values  of  the  thickness  in  meters  along 
the  x-axis  and  the  diagonal. 

Clamped  sandwich  panels  subjected  to  a  point  load 

Weight  saving  for  simply  supported  sandwiches  subjected  to  a  point  load, 
is  shown  in  Figure  106.  The  weight  saving  is  a  function  of  the  aspect  ratio, 
the  maximum  value  being  96  %  for  b/a=10.  The  optimum  sublaminate  also 
varies  in  the  function  of  the  aspect  ratio: 

-  [0/90]  for  b/a  between  1  and  1.5, 

.  [O5/9O2]  for  b/a  between  1.5  and  5, 

.  [0]  for  b/a  higher  than  5. 

The  quasi-isotropic  sublaminate  [0/45/-45/90]  presents  low-weight  saving 
with  respect  to  the  optimum  sublaminate,  the  maximum  value  being  94  %  for 
b/a=2. 

Since  no  general  formula  for  expressing  the  optimum  laminate  thickness 
has  been  found,  the  following  particular  case  has  been  calculated: 

.  a=7.87  in  (0.2  m), 
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.  P=1.125  E4  lb  (5  E4  N). 


Figures  107  through  109  show  the  values  of  the  thickness  in  meters  along 
the  x-axis  and  the  diagonal. 
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4.  OPTIMIZATION  OF  LAMINATED  COMPOSITE  PLATES 
SUBJECTED  TO  BUCKLING 


This  section  analyzes  the  case  dealing  with  buckling  loads.  The  use  of 
laminated  composite  materials  in  aircrafts  and  space  structures  offers  a 
unique  opportunity  for  obtaining  weight  savings  over  conventional 
constructions.  Also,  composite  materials  adds  versatility  to  the  design  process 
by  allowing  the  structure  to  be  tailored  to  meet  the  design  criteria.  One  such 
design  criterion  is  the  prevention  of  compressive  and  shear  buckling  in 
laminated  composite  plates. 

First,  the  analytical  formulation  of  buckling  of  plates  is  reported.  As  usual, 
the  theory  of  plates  requires  a  numerical  technique  in  order  to  simulate 
different  boundary  conditions  and  types  of  loading.  The  same  finite  element 
used  in  Section  3  (Optimization  of  laminated  composite  plates  subjected  to 
transverse  loads)  has  been  applied  in  this  work. 

Finally,  results  of  one-  and  two-dimensional  laminated  composite  plates  have 
been  described.  Optima  angles  of  orientation  of  the  fiber  and  critical  uniform 
compressive  and  shear  loads  are  given  for  different  boundary  conditions  and 
types  of  loading.  Relative  weights  with  respect  to  aluminium  and  the  quasi¬ 
isotropic  configuration  are  also  reported. 


4.1  ANALYTICAL  FORMULATION 


The  strain  energy  of  a  single  plate  element  subject  to  in-plane  stresses  can 
be  written  as 

V=VB+VS 

(4.1) 


62 


where  Vs  is  the  strain  energy  owing  to  in-plane  stresses  acting  on  the  second 
order  strains,  and  VB  is  the  strain  energy  owing  to  plate  bending.  Vs  can  be 
expressed  as 


te  (~)+Oy  (np)2+  2txy 
dx  dy 


dw  dw 
dx  dy 


]  t  dx  dy 


or  written  in  matrix  form  as 

dw 

5"  ^llfltdxdy 

dx  dy  ^xy  dw 

¥ 


(4.2) 


(4.3) 


Now,  the  lateral  displacement  w  can  be  written  as 
w  =  [N]  {8e} 


(4.4) 


where  N  is  the  shape  function  vector  and  {8e}  is  a  vector  of  nodal 
deflections.  Writing 


dw 


dN 


(  dx  ]  =  [  _  ]  (S')  =  [SJ  (S') 
2*  3N 


dy 


dy 


(4.5) 


the  following  expression  can  be  obtained 


Vs  =  -1 
2 


mT[S]T[c]  [S]  {5e}  tdx  dy 


(4.6) 
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On  differentiating  with  respect  to  the  nodal  displacements  this  component 
of  strain  energy  gives  a  "geometric  stiffness" 

9V5 

3Si 

[  .  ]  =  -[ks]e{8e} 


av6 

asN 

with 


[ks]e 


[S]T[a]  [Sjtdxdy 


(4.7) 


(4.8) 


For  the  whole  plate,  the  element  matrices  can  be  assembled: 

( [  KB]  -  [  Ks]  )  {5}  =  {Q}  (4.9) 

If  all  the  in-plane  stresses  can  be  increased  by  a  factor  A.,  this  will  increase 
the  geometric  stiffness  proportionately.  It  is  possible,  therefore,  to  find  a  A 
which,  with  no  external  loading  Q,  gives 

( [  Kb]  -  A  [  Ks]  )  {8}  =  0  (4.10) 

that  expresses  an  eigenvalue  problem. 


4.2  DESCRIPTION  OF  RESULTS 
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In  this  part  of  the  report,  some  results  of  laminates  composite  plates 
subjected  to  compression,  shear  and  combined  loads  are  presented.  Only  thin 
plates  have  been  analyzed. 

In  the  case  of  one-dimensional  plates,  all  the  cases  reported  are  for  an 
angle  of  orientation  of  the  fibers  of  0°.  Critical  loads,  weight  savings  for  flat 
plates  and  different  cross  sections  are  presented. 

In  the  chapter  about  two-dimensional  plates,  the  study  of  variable  thickness 
plates  has  not  generated  interesting  results.  The  only  author  that  has  published 
in  this  area  was  Ashton  [75],  who  obtained  results  for  linearly  tapered  plates 
subjected  to  uniaxial  loading,  but  no  weight  saving  was  reported.  Instead  of 
that,  very  interesting  results  have  appeared  when  the  angle  of  orientation  of 
the  fiber  varies. 


4.2.1  Model  assumptions  and  method  of  analysis 


The  model  used  presents  the  following  assumptions: 

The  analysis  is  static. 

The  theory  applied  is  valid  for  small  deformations. 

The  material  used  is  supposed  to  be  elastic. 

The  finite  element  used  here  is  based  on  the  higher  order  shear 
theory [62]  and  the  penalty  function  method.  A  2-D  model  has  been 
applied  for  both  one-  and  two-dimensional  laminated  composite  plates. 

The  design  criterion  used  is  the  critical  compressive  force  to  prevent 
the  plate  from  buckling. 
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4.2.2  One-dimensional  laminated  composite  plates 
The  following  points  analyze  the  problem: 

.  The  objective  is  to  get  the  minimum  weight  structure  by  using  critical 
compressive  force  as  a  design  criterion. 

.  The  laminate  used  is  [0]n, because  it  is  the  optimum  for  one-dimensional 
laminated  composite  plates  subjected  to  uniform  uniaxial  loads. 

.  Boundary  conditions  are  restricted  to  : 

-  simply  supported  plates, 

-  clamped  plates, 

-  cantilever  plates. 

.  Uniform  uniaxial  compression  loading  has  been  applied. 

.  Plate  aspect  ratios  between  1  and  10  have  been  considered. 

.  The  mesh  was  composed  by  400  nodes. 

.  The  material  used  was  T300/N5208. 

Figure  1 10  shows  a  definition  of  coordinates  axis  for  one-dimensional  plates 
subjected  to  uniform  uniaxial  compression  loads. 

One-dimensional  simply  supported  plates  subjected  to  a  uniaxial  compression 
load 

Figures  111  and  112  show  critical  loads  for  one-dimensional  simply 
supported  plates  subjected  to  uniform  uniaxial  compression  loads.  Figure  1 1 1 
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also  shows  the  optimum  angle.  Since  this  case  is  one-dimensional,  this  angle 
is  0°.  A  comparison  between  the  optimum  angle  and  the  quasi -isotropic 
laminate  is  shown  in  Figure  112.  This  laminate  presents  very  low  critical 
loads  for  low-aspect  ratios. 

Optima  configurations  for  one-dimensional  simply  supported  plates 
subjected  to  uniform  uniaxial  compression  loads  are  depicted  in  Figures  113 
and  114.  For  a  maximum  height  of  twice  the  thickness,  the  triangular 
configuration  presents  a  higher  weight  saving  than  the  other  two  (Figure 
1 13).  The  T  configuration  is  the  lightest  if  the  height  is  ten  times  the  thickness 
(Figure  114). 

One- dimensional  clamped  plates  subjected  to  a  uniaxial  compression  load 

Figures  115  through  118  refer  to  one-dimensional  clamped  plates  subjected 
to  uniform  uniaxial  compression  loads.  A  comparison  between  the  optimum 
angle,  0°  (Figure  115)  and  the  quasi-isotropic  laminate  is  shown  in  Figure 
1 16.  The  critical  load  for  this  case  is  four  times  higher  than  it  was  in  the  case 
of  simply  supported  plates.  The  reason  for  this  is  that  a  clamped  plate  is 
stiffer  than  a  simply  supported  one  and  thus,  the  capability  of  supporting 
uniform  uniaxial  compression  loads  is  also  higher.  In  the  chapter  of  optima 
configurations,  similar  results  to  the  case  of  simply  supported  plates  are 
reported  (Figures  117  and  118). 

One-dimensional  cantilever  plates  subjected  to  a  uniaxial  compression  load 

Finally,  one-dimensional  cantilever  plates  subjected  to  uniform  uniaxial 
compression  loads  are  analyzed  in  Figures  119  through  122.  The  critical  load 
for  this  case  is  4  times  lower  than  it  was  in  the  case  of  simply  supported 
plates,  and  16  times  lower  than  in  the  case  of  clamped  plates  (Figure  119).  A 
comparison  between  the  optimum  angle  and  the  quasi-isotropic  laminate  is 
shown  in  Figure  1 20.  In  the  chapter  of  optima  configurations,  similar  results 
to  the  two  boundary  conditions  treated  above  are  reported  (Figures  121  and 
122). 
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4.2.3  Two-dimensional  laminated  composite  plates 


The  following  points  analyze  the  problem: 

.  The  objective  is  to  get  the  minimum  weight  structure  by  using  critical 
compressive  force  as  a  design  criterion. 

.  Only  symmetric  and  balanced  laminates  have  been  considered. 

.  Boundary  conditions  are  restricted  to  : 

-  simply  supported  plates  along  the  four  edges  of  the  plate,  and 

-  clamped  plates  along  the  four  edges  of  the  plate. 

.  Six  types  of  loading  have  been  applied: 

-  uniform  uniaxial  compression  load, 

-  uniform  biaxial  compression  load  (Ny=Nx/2), 

-  uniform  biaxial  compression  load  (Ny=Nx), 

-  uniform  biaxial  compression  load  (Ny=2*Nx), 

-  uniform  shear  load, 

-  combined  load:  uniform  uniaxial  and  uniform  shear  load. 

.  The  mesh  was  composed  by  400  nodes. 

.  The  material  used  is  T300/N5208. 

Figure  123  shows  a  definition  of  coordinates  axis  for  two-dimensional 
plates  subjected  to  uniform  uniaxial,  biaxial,  shear  and  combined  compression 
loads. 

Two-dimensional  simply  supported  plates  subjected  to  a  uniaxial  compression 
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In  Figures  124  and  125,  critical  loads  are  reported  for  two-dimensional 
simply  supported  plates  subjected  to  uniform  uniaxial  compression  loads. 
Figure  124  also  shows  the  optimum  angle.  Both  critical  load  and  optimum 
angle  are  a  function  of  the  plate  aspect  ratio,  though  they  present  constant 
values  for  aspect  ratios  higher  than  four.  A  comparison  between  the  optimum 
configuration  and  the  quasi-isotropic  laminate  is  shown  in  Figure  125. 

A  comparison  between  critical  loads  for  angle-ply  and  quasi-isotropic 
laminate  in  function  of  the  angle  of  fiber  orientation  is  carried  out  in  Figure 
126.  For  angles  between  27  and  66  degrees,  the  critical  load  for  angle-ply 
laminates  is  higher  than  for  quasi-isotropic  laminates. 

Figure  127  shows  normalized  weights  of  aluminum  and  quasi-isotropic 
T300/N5208  to  the  optimum  laminate  of  T300/N5208. 


Two-dimensional  simply  supported  plates  subjected  to  a  biaxial  compression 
(Ny=Nj2) 

Figures  128  and  129  report  critical  loads  for  two-dimensional  simply 
supported  plates  subjected  to  uniform  biaxial  compression  loads  (Ny=Nx/2). 

Figure  128  also  shows  the  optimum  angle.  Both  critical  load  and  optimum 
angle  are  a  function  of  the  plate  aspect  ratio,  though  they  present  constant 
values  for  aspect  ratios  higher  than  two.  A  comparison  between  the  optimum 
configuration  and  the  quasi-isotropic  laminate  is  shown  in  Figure  129. 

A  comparison  between  critical  loads  for  angle-ply  and  quasi-isotropic 
laminate  in  function  of  the  angle  of  fiber  orientation  is  carried  out  in  Figure 
130.  For  angles  between  27  and  66  degrees,  the  critical  load  for  angle-ply 
laminates  is  higher  than  for  quasi-isotropic  laminates. 

Figure  131  shows  normalized  weights  of  aluminum  and  quasi-isotropic 
T300/N5208  to  the  optimum  laminate  of  T300/N5208. 


69 


Two-dimensional  simply  supported  plates  subjected  to  a  biaxial  compression 
(Ny  =  Nx) 

Figures  132  and  133  reports  critical  loads  for  two-dimensional  simply 
supported  plates  subjected  to  uniform  biaxial  compression  loads  (Ny=Nx). 

Figure  132  also  shows  the  optimum  angle.  Both  critical  load  and  optimum 
angle  are  a  function  of  the  plate  aspect  ratio,  though  they  present  constant 
values  for  aspect  ratios  higher  than  three.  A  comparison  between  the 
optimum  configuration  and  the  quasi-isotropic  laminate  is  shown  in  Figure 
133. 


A  comparison  between  critical  loads  for  angle-ply  and  quasi-isotropic 
laminate  in  function  of  the  angle  of  fiber  orientation  is  carried  out  in  Figure 
134.  For  angles  between  27  and  66  degrees,  the  critical  load  for  angle-ply 
laminates  is  higher  than  for  quasi-isotropic  laminates. 

Figure  135  reports  normalized  weights  of  aluminum  and  quasi-isotropic 
T300/N5208  to  the  optimum  laminate  of  T300/N52 08. 

Two-dimensional  simply  supported  plates  subjected  to  a  biaxial  compression 
(Ny-2*NX) 

Figures  136  and  137  show  the  critical  loads  for  two-dimensional  simply 
supported  plates  subjected  to  uniform  biaxial  compression  loads  (Ny=2*Nx). 

Figure  136  also  shows  the  optimum  angle.  Both  critical  load  and  optimum 
angle  are  a  function  of  the  plate  aspect  ratio,  though  they  present  constant 
values  for  aspect  ratios  higher  than  two.  A  comparison  between  the  optimum 
configuration  and  the  quasi-isotropic  laminate  is  shown  in  Figure  137. 

A  comparison  between  critical  loads  for  angle-ply  and  quasi-isotropic 
laminate  in  function  of  the  angle  of  fiber  orientation  is  carried  out  in  Figure 
138.  For  angles  between  27  and  66  degrees,  the  critical  load  for  angle-ply 
laminates  is  higher  than  for  quasi-isotropic  laminates. 
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Figure  139  shows  normalized  weights  of  aluminum  and  quasi-isotropic 
T300/N5208  to  the  optimum  laminate  of  T300/N5208. 

Two-dimensional  clamped  plates  subjected  to  a  uniaxial  compression 

Figures  140  and  141  report  critical  loads  for  two-dimensional  clamped 
plates  subjected  to  uniform  uniaxial  compression  loads.  Figure  140  also 
shows  the  optimum  angle.  Both  critical  load  and  optimum  angle  are  a 
function  of  the  plate  aspect  ratio.  A  comparison  between  the  optimum 
configuration  and  the  quasi-isotropic  laminate  is  shown  in  Figure  141. 

A  comparison  between  critical  loads  for  angle-ply  and  quasi-isotropic 
laminate  in  function  of  the  angle  of  fiber  orientation  is  carried  out  in  Figure 
142.  For  angles  lower  than  42  degrees,  the  critical  load  for  angle-ply 
laminates  is  higher  than  for  quasi-isotropic  laminates. 

Figure  143  reports  normalized  weights  of  aluminum  and  quasi-isotropic 
T300/N5208  to  the  optimum  laminate  of  T300/N5208. 

Two-dimensional  clamped  plates  subjected  to  a  biaxial  compression 
(Ny=Nx/2) 

Figures  144  and  145  reports  critical  loads  for  two-dimensional  clamped 
plates  subjected  to  uniform  biaxial  compression  loads  (Ny=Nx/2).  Figure  144 

also  shows  the  optimum  angle.  Both  critical  load  and  optimum  angle  are  a 
function  of  the  plate  aspect  ratio.  A  comparison  between  the  optimum 
configuration  and  the  quasi-isotropic  laminate  is  shown  in  Figure  145. 

Comparison  between  critical  loads  for  angle-ply  and  quasi-isotropic 
laminate  in  function  of  the  angle  of  fiber  orientation  is  carried  out  in  Figure 
146.  For  angles  between  21  and  40  degrees,  the  critical  load  for  angle -ply 
laminates  is  higher  than  for  quasi-isotropic  laminates. 

Figure  147  reports  normalized  weights  of  aluminum  and  quasi-isotropic 
T300/N5208  to  the  optimum  laminate  of  T300/N5208. 
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Two-dimensional  clamped  plates  subjected  to  a  biaxial  compression  (Ny=Nx) 


Figures  148  and  149  report  critical  loads  for  two-dimensional  clamped 
plates  subjected  to  uniform  biaxial  compression  loads  (Ny=Nx).  Figure  148 

also  shows  the  optimum  angle.  Both  critical  load  and  optimum  angle  are  a 
function  of  the  plate  aspect  ratio.  A  comparison  between  the  optimum 
configuration  and  the  quasi-isotropic  laminate  is  shown  in  Figure  149. 

Comparison  between  critical  loads  for  angle-ply  and  quasi-isotropic 
laminate  in  function  of  the  angle  of  fiber  orientation  is  carried  out  in  Figure 
150.  For  angles  between  43  and  47  degrees,  the  critical  load  for  angle-ply 
laminates  is  higher  than  for  quasi-isotropic  laminates. 

Figure  151  shows  normalized  weights  of  aluminum  and  quasi-isotropic 
T300/N5208  to  the  optimum  laminate  of  T300/N5208. 


Two-dimensional  clamped  plates  subjected  to  a  biaxial  compression 
(Ny=2*Nx) 

Figures  152  and  153  report  critical  loads  for  two-dimensional  clamped 
plates  subjected  to  uniform  biaxial  compression  loads  (Ny=2*Nx).  Figure  152 

also  .shows  the  optimum  angle.  Both  critical  load  and  optimum  angle  are  a 
function  of  the  plate  aspect  ratio.  A  comparison  between  the  optimum 
configuration  and  the  quasi-isotropic  laminate  is  shown  in  Figure  153. 

A  comparison  between  critical  loads  for  angle-ply  and  quasi-isotropic 
laminate  in  function  of  the  angle  of  fiber  orientation  is  carried  out  in  Figure 
154.  For  angles  between  52  and  66  degrees,  the  critical  load  for  angle-ply 
laminates  is  higher  than  for  quasi-isotropic  laminates. 

Figure  155  shows  normalized  weights  of  aluminum  and  quasi-isotropic 
T300/N5208  to  the  optimum  laminate  of  T300/N5208. 
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Two-dimensional  simply  supported  plates  subjected  to  a  shear  load. 

Figures  156  and  157  reports  critical  loads  for  two-dimensional  simply 
supported  plates  subjected  to  uniform  shear  loads.  Figure  156  also  shows  the 
optimum  angle.  Both  critical  load  and  optimum  angle  are  a  function  of  the 
plate  aspect  ratio.  A  comparison  between  the  optimum  configuration  and  the 
quasi-isotropic  laminate  is  shown  in  Figure  157. 

A  comparison  between  critical  loads  for  angle-ply  and  quasi-isotropic 
laminate  in  function  of  the  angle  of  fiber  orientation  is  carried  out  in  Figure 
158.  For  angles  between  24  and  66  degrees,  the  critical  load  for  angle-ply 
laminates  is  higher  than  for  quasi-isotropic  laminates. 

Figure  159  reports  normalized  weights  of  aluminum  and  quasi-isotropic 
T300/N5208  to  the  optimum  laminate  of  T300/N5208. 


Two-dimensional  clamped  plates  subjected  to  a  shear  load. 

Figures  160  and  161  report  critical  loads  for  two-dimensional  clamped 
plates  subjected  to  uniform  shear  loads.  Figure  160  also  shows  the  optimum 
angle.  Both  critical  load  and  optimum  angle  are  a  function  of  the  plate  aspect 
ratio.  A  comparison  between  the  optimum  configuration  and  the  quasi¬ 
isotropic  laminate  is  shown  in  Figure  161. 

A  comparison  between  critical  loads  for  angle-ply  and  quasi-isotropic 
laminate  in  function  of  the  angle  of  fiber  orientation  is  carried  out  in  Figure 
162.  For  angles  between  26  and  64  degrees,  the  critical  load  for  angle-ply 
laminates  is  higher  than  for  quasi-isotropic  laminates. 

Figure  163  shows  normalized  weights  of  aluminum  and  quasi-isotropic 
T300/N5208  to  the  optimum  laminate  of  T300/N5208. 

Two-dimensional  simply  supported  and  clamped  plates  subjected  to  a 
combined  load  (uniaxial  compression  and  shear) 
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Finally,  Figures  164  and  165  represent  the  buckling  parameters  for  simply 
supported  and  clamped  plates  subjected  to  combined  uniform  uniaxial 
compression  and  shear,  respectively. 
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5.  CONCLUDING  REMARKS 


A  study  on  optimization  of  laminated  composite  plates  has  been  presented. 
Plates  subjected  to  transverse  loads  and  buckling  have  been  treated.  Optima 
configurations  for  different  boundary  conditions  and  types  of  loading  have 
been  reported. 

In  the  chapter  of  transverse  loads,  the  following  conclusions  can  be  drawn: 

.  The  variable  thickness  effect  has  a  remarkable  influence  on  the  the 
behavior  of  a  tapered  laminated  composite  plate.  The  distributions  of 
through-thickness  stresses  are  very  sensitive  to  the  thickness  ratio  (t/tl) 
and  the  angle  of  variation  of  thickness.  Especially,  the  interlaminar  shear 
stress  reaches  very  high  values  in  the  areas  of  change  of  thickness.  The 
peak  value  can  be  controlled  by  modifying  some  parameters:  thickness 
ratio,  angle  of  variation  of  thickness,  type  of  loading....  The  distributions 
of  interlaminar  stresses  vary  in  function  of  the  sign  of  the  bending 
moment. 

.  The  failure  mechanisms  reported  in  variable-thickness  composite  plates 
can  be  grouped  in  three  general  modes: 

1)  Interlaminar  shear  mode:  It  is  due  to  interlaminar  shear  stresses,  and 
appears  in  the  tapered  surfaces  when  the  change  of  thickness  is  sharp 
(high  angles  of  variation  of  thickness) .  It  can  also  appear  in  the 
midplane,  in  thick  plates. 

2)  Bending  modes:  They  are  due  to  bending  stresses.  These  modes  are 
usual  in  thin  plates,  very  low  angles  of  variation  of  thickness  and  low 
thickness  ratios.  They  are  the  typical  failure  modes  in  thin  untapered 
plates. 
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3)  Mixed  modes:  A  compression-interlaminar  normal  mode  has  been 
found  in  the  area  of  application  of  the  load  in  the  fatigue  testing.  This 
mode  does  not  appear  in  the  4-point  bending  test. 

.  For  uniformly  distributed  loads,  optimal  configurations  have  been  shown. 
Some  general  formulas  have  been  presented  for  designing  variable¬ 
thickness  laminated  composite  plates  (3.47  )  to  (3.52).  High  weight  savings 
have  also  been  reported. 

.  For  point  loads,  we  have  found  a  general  formula  only  for  one¬ 
dimensional  plates.  Hence,  a  particular  case  has  been  analyzed  for  two- 
dimensional  plates.  For  both  cases,  one-  and  two-dimensional  plates,  the 
optimal  configurations  have  been  reported.  Very  high  weight  saving  have 
also  been  presented. 

From  the  chapter  of  buckling  loads,  the  following  conclusions  can  be 

drawn: 

.  Optimal  configurations  are  presented  for  one-dimensional  plates.  Critical 
compressive  forces  are  also  represented.  Different  variable-thickness 
solutions  are  depicted,  and  their  weight  savings  given  (from  15  to  66%). 

.  Finally,  optimal  two-dimensional  plates  subjected  to  uniform  uni-,biaxial 
compression,  shear,  and  combined  loads  are  shown.  Optima  angles  of 
orientation  of  fibers,  critical  compressive  and  shear  forces  have  been 
presented  for  all  cases,  and  relative  weights  with  respect  to  aluminium 
given.  The  highest  weight  saving  have  been  reported  for  low  aspect 
ratios. 
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Figure  1.  Schematic  representation  of  a  variable-thickness  composite  plate. 


Figure  2.  Representation  of  the  optimization  procedure. 
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Figure  1 1 .  Distribution  of  05  for  a  two-sides  tapered  plate  and  dM/3x  <  0. 


Figure  12.  Distribution  of  G\  for  a  one-side  tapered  plate  and  3M/8x  <  0. 
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Figure  15.  Distribution  of  Oj  for  a  one-side  tapered  plate,  ti/t=0.3,  M/Pt=5  and 
9M/dx  >  0. 
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Figure  16.  Distribution  of  03  for  a  one-side  tapered  plate,  ti/t=0.3,  M/Pt=5  and 
dM/dx  >0. 
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Figure  17.  Distribution  of  <75  for  a  one-side  tapered  plate,  ti/t=0.3,  M/Pt=5  and 
dM/dx  >  0. 


Figure  18.  Distribution  of  Oj  for  a  one-side  tapered  plate,  tj/t=0.2,  M/Pt=5  and 
dM/dx  >  0. 
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Figure  21.  Distribution  of  oi  for  a  one-side  tapered  plate,  t]/t=0.1,  M/Pt=5  and 
8M/8x  >0. 
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Figure  22.  Distribution  of  03  for  a  one-side  tapered  plate,  ti/t=0.1,  M/Pt=5  and 
8M/8x  >0. 
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Figure  23.  Distribution  of  05  for  a  one-side  tapered  plate,  tj/t=0.1,  M/Pt=5  and 
dM/dx  >0. 


Figure  24.  ibree-point  bending  test. 


Figure  25.  Delamination  failure  in  the  midplane. 


Figure  26.  Delamination  failure  in  the  tapered  surface 
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'  21 '  SEM  Ph0t0graPh  of  delamination  fracture  surface. 
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Figure  29.  Photograph  showing  a  mixed  compression-interlaminar  normal  failu. 
mode. 
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Figure  30.  SEM  photograph  of  mixed  mode  fracture  surface 
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Figure  31.  SEM  image  of  top  of  mixed  failure  mode  surface.  Crushing  due  to  a 
compressive  interlaminar  stress  is  detected. 


BEAM  ‘  - 

0.008  1 0KU  XI.  000  1 0 k- m  W D 3 1 


Figure  32.  SEM  photograph  of  the  upper  part  of  the  mixed  failure  mode  surface. 
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Fieure  33.  SbM  photograph  of  the  P 
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Figure  34.  SEM  image  of  die  inner  part 


of  the  mixed  failure  mode  surface. 


Figure  35.  SEM  image  of  the  bottom  of  the  mixed  failure  mode  surface. 


Figure  36.  Photograph  showing  specimen  OLCP_7001. 
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Figure  39.  Image  of  failure  in  specimen  OLCP_7001. 


Figure  40.  Representation  of  specimen  OLCPJ7001  and  position  of  gauges. 
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Figure  41.  Photograph  showing  specimen  OLCP_7002. 
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Figure  42.  Measurements  of  specimen  OLCP_7002  in  inches  (meters). 


Figure  43.  Image  of  failure  in  specimen  OLCP_7002. 


Figure  44.  Representation  of  specimen  OLCP_7002  and  position  of  gauges. 
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Figure  45.  Photograph  showing  specimen  OLCP_7003. 
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Figure  46.  Measurements  of  specimen  OLCPJ7003  in  inches  (meters). 


Figure  47.  Image  of  failure  in  specimen  OLCP_7003. 
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Figure  48.  Representation  of  specimen  OLCP_7003  and  position  of  gauges. 


Figure  49.  Representation  of  sections  A  A  and  BB. 
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Figure  50.  Distribution  of  £5  through  the  thickness  in  section  AA 
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Figure  51.  Distribution  of  £5  through  the  thickness  in  section  BB 


Figure  52.  Description  of  structure  used  to  compare  both  2-D  plane  strain  model 
and  shear  deformation  plate  model. 
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Figure  53.  Distribution  of  Oi  through  the  thickness  from  a  2-D  plane  strain  model 
and  a  shear  deformation  plate  model. 
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Figure  54.  Distribution  of  05  through  the  thickness  from  a  2-D  plane  strain  model 
and  a  shear  deformation  plate  model. 
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Figure  55.  Variation  of  interlaminar  normal  modulus  in  function  of  m  and  n  in  a 
laminate  T300/N5208  [90m  /0n]  s. 
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Figure  56.  Variation  of  interlaminar  shear  modulus  in  function  of  (J>  in  a  laminate 
T300/N5208  [+<}>  /-(j)  ]  s. 


Figure  57.  Definition  of  coordinates  axis  and  key  to  colors  for  one-dimensional 
plates. 
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Figure  59.  Representation  of  optimal  thick  one-dimensional  plates  subjected  to  a 
uniform  distributed  transverse  load. 
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Figure  60.  Representation  of  optimal  one-dimensional  sandwich  panels  subjected 
to  a  uniform  distributed  transverse  load. 
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Figure  62.  Representation  of  optimal  thick  one-dimensional  plates  subjected  to  a 
point  transverse  load. 
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Figure  63.  Representation  of  optimal  one-dimensional  sandwich  panels  subjected  to 
a  point  transverse  load. 
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Figure  64.  Definition  of  coordinate  axis  and  key  to  colors  for  two-dimensional 
plates  subjected  to  a  transverse  load. 
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Figure  65.  Weight  saving  and  normalized  deflection  for  a  thin,  simply  supported 
plate  subjected  to  a  uniform  distributed  transverse  load. 
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Figure  66.  Representation  of  an  optimal  thin,  simply  supported  plate  subjected  to  a 
uniform  distributed  transverse  load  for  aspect  ratios  between  1  and  1.75. 
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Figure  67.  Representation  of  an  optimal  thin,  simply  supported  plate  subjected  to  a 
uniform  distributed  transverse  load  for  aspect  ratios  between  1.75  and  2.25. 
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Figure  68.  Representation  of  an  optimal  thin,  simply  supported  plate  subjected  to  a 
uniform  distributed  transverse  load  for  aspect  ratios  between  2.25  and  3.5. 
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Figure  69.  Representation  of  an  optimal  thin,  simply  supported  plate  subjected  to  a 
uniform  distributed  transverse  load  for  aspect  ratios  higher  than  3.5. 
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Figure  70.Weight  saving  and  normalized  deflection  for  a  thin,  clamped  plate 
subjected  to  a  uniform  distributed  transverse  load. 
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Figure  7 1 .  Representation  of  an  optimal  thin,  clamped  plate  subjected  to  a  uniform 
distributed  transverse  load  for  aspect  ratios  between  1  and  1 .5. 
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Figure  72.  Representation  of  an  optimal  thin,  clamped  plate  subjected  to  a  uniform 
distributed  transverse  load  for  aspect  ratios  between  1.5  and  4.4. 
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Figure  73.  Representation  of  an  optimal  thin,  clamped  plate  subjected  to  a  uniform 
distributed  transverse  load  for  aspect  ratios  higher  than  4.4. 
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Figure  74.  Weight  saving  and  normalized  deflection  for  a  thick,  simply  supported 
plate  subjected  to  a  uniform  distributed  transverse  load. 
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Figure  75.  Representation  of  an  optimal  thick,  simply  supported  plate  subjected  to 
a  uniform  distributed  transverse  load  plates  for  aspect  ratios  between  1  and  3.5. 
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Figure  76.  Representation  of  an  optimal  thick,  simply  supported  plate  subjected  to 
a  uniform  distributed  transverse  load  for  aspect  ratios  higher  than  3.5. 
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Figure  77.  Weight  saving  and  normalized  deflection  for  a  thick,  clamped  plate 
subjected  to  a  uniform  distributed  transverse  load. 
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Figure  78.  Representation  of  an  optimal  thick,  clamped  plate  subjected  to  a 
uniform  distributed  transverse  load  for  aspect  ratios  between  1  and  3.5. 
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Figure  79.  Representation  of  an  optimal  thick,  clamped  plate  subjected  to  a 
uniform  distributed  transverse  load  for  aspect  ratios  higher  than  3.5. 
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Figure  80.  Weight  saving  and  normalized  deflection  for  a  simply  supported 
sandwich  panel  subjected  to  a  uniform  distributed  transverse  load. 


Figure  81.  Representation  of  an  optimal  simply  supported  sandwich  panel  subjected 
to  a  uniform  distributed  transverse  load  for  aspect  ratios  between  1  and  1.5. 
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Figure  82.  Representation  of  an  optimal  simply  supported  sandwich  panel  subjected 
to  a  uniform  distributed  transverse  load,  for  aspect  ratios  higher  than  1.5. 
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Figure  83.  Weight  saving  and  normalized  deflection  for  a  clamped  sandwich  panel 
subjected  to  a  uniform  distributed  transverse  load. 
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Figure  84.  Representation  of  an  optimal  clamped  sandwich  panel  subjected  to  a 
uniform  distributed  transverse  load  for  aspect  ratios  between  1  and  1.5. 


Figure  85.  Representation  of  an  optimal  clamped  sandwich  panel  subjected  to  a 
uniform  distributed  transverse  load  for  aspect  ratios  between  1.5  and  2.5. 
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Figure  86.  Representation  of  an  optimal  clamped  sandwich  panel  subjected  to  a 
uniform  distributed  transverse  load  for  aspect  ratios  higher  than  2.5. 
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Figure  87.  Weight  saving  for  a  thin,  simply  supported  plate  subjected  to  a  point 
transverse  load. 
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Figure  88.  Representation  of  an  optimal  thin,  simply  supported  plate  subjected  to  a 
point  transverse  load  for  aspect  ratios  between  1  and  2. 


Figure  89.  Representation  of  an  optimal  thin,  simply  supported  plate  subjected  to  a 
point  transverse  load  for  aspect  ratios  between  2  and  5. 
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Figure  90.  Representation  of  an  optimal  thin,  simply  supported  plate  subjected  to  a 
point  transverse  load  for  aspect  ratios  higher  than  5. 
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Figure  91.  Weight  saving  for  a  thin,  clamped  plate  subjected  to  a  point  transverse 


Figure  92.  Representation  of  an  optimal  thin,  clamped  plate  subjected  to  a  point 
transverse  load  for  aspect  ratios  between  1  and  1.5. 


Figure  93.  Representation  of  an  optimal  thin,  clamped  plate  subjected  to  a  point 
transverse  load  for  aspect  ratios  between  1.5  and  5. 
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Figure  94.  Representation  of  an  optimal  thin,  clamped  plate  subjected  to  a  point 
transverse  load  for  aspect  ratios  higher  than  5. 
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Figure  95.  Weight  saving  for  a  thick,  simply  supported  plate  subjected  to  a  point 
transverse  load. 


132 


b 

1  <  -  <  2  25 

”  u 

OPTIMUM  SUBLAMINATE 
I45/-45] 

a  =  0.01  m 
P  =4.0  E3  N 


Figure  96.  Representation  of  an  optimal  thick,  simply  supported  plate  subjected  to 
a  point  transverse  load  for  aspect  ratios  between  1  and  2.25. 
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Figure  97.  Representation  of  an  optimal  thick,  simply  supported  plate  subjected  to 
a  point  transverse  load  for  aspect  ratios  higher  than  2.25. 
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Figure  98.  Weight  saving  for  a  thick,  clamped  plate  subjected  to  a  point 
transverse  load. 


Figure  99.  Representation  of  an  optimal  thick,  clamped  plate  subjected  to  a  point 
transverse  load  for  aspect  ratios  between  1  and  2.5. 


Figure  100.  Representation  of  an  optimal  thick,  clamped  plate  subjected  to  a  point 
transverse  load  for  aspect  ratios  between  2.5  and  5. 
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Figure  101.  Representation  of  an  optimal  thick,  clamped  plate  subjected  to  a  point 
transverse  load  for  aspect  ratios  higher  than  5. 
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Figure  102.  Weight  saving  for  a  simply  supported  sandwich  panel  subjected  to  a 
point  transverse  load. 


Figure  103.  Representation  of  an  optimal  simply  supported  sandwich  panel 
subjected  to  a  point  transverse  load  for  aspect  ratios  between  1  and  2. 
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Figure  104.  Representation  of  an  optimal  simply  supported  sandwich  panel 
subjected  to  a  point  transverse  load  for  aspect  ratios  between  2  and  5. 
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Figure  105.  Representation  of  an  optimal  simply  supported  sandwich  panel 
subjected  to  a  point  transverse  load  for  aspect  ratios  higher  than  5. 
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Figure  106.  Weight  saving  for  a  clamped  sandwich  panel  subjected  to  a  point 
transverse  load. 


Figure  107.  Representation  of  an  optimal  clamped  sandwich  panel  subjected  to  a 
point  transverse  load  for  aspect  ratios  between  1  and  1.5. 


138 


b 

«5<  5 

OPTIMUM  SUBLAMINATE 
[05/902] 


0:1m 
c  =  0.1  m 
P  =  2.0  E4N 


0  2  x  1 


1 


Figure  108.  Representation  of  an  optimal  clamped  sandwich  panel  subjected  to  a 
point  transverse  load  for  aspect  ratios  between  1.5  and  5. 
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Figure  109.  Representation  of  an  optimal  clamped  sandwich  panel  subjected  to  a 
point  transverse  load  for  aspect  ratios  higher  than  5. 


Figure  110.  Definition  of  coordinate  axis  for  one-dimensional  plates  subjected  to  a 
uniform  uniaxial  compression  load. 
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Figure  111.  Critical  load  and  optimum  angle  for  an  one-dimensional  simply 
supported  plate. 
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Figure  112.  Comparison  between  critical  loads  of  optimum  configuration  and 
quasi-isotropic  laminate  (one-dimensional  simply  supported  plates). 
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Figure  1 13.  Weight  saving  for  different  cross  sections  whose  height  is  twice  the 
laminate  thickness  (one-dimensional  simply  supported  plates). 
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Figure  114.  Weight  saving  for  different  cross  sections  whose  height  is  ten  times 
the  laminate  thickness  (one-dimensional  simply  supported  plates). 


Figure  115.  Critical  load  and  optimum  angle  for  an  one-dimensional  clamped 
plate. 
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Figure  116.  Comparison  between  critical  loads  of  optimum  configuration  and 
quasi-isotropic  laminate  (one-dimensional  clamped  plates). 
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Figure  117.  Weight  saving  for  different  cross  sections  whose  height  is  twice  the 
laminate  thickness  (one-dimensional  clamped  plates). 
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Figure  118.  Weight  saving  for  different  cross  sections  whose  height  is  ten  times 
the  laminate  thickness  (one-dimensional  simply  supported  plates). 


Figure  119.  Critical  load  and  optimum  angle  for  an  one-dimensional  cantilever 
plate. 
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Figure  120.  Comparison  between  critical  loads  of  optimum  configuration  and 
quasi-isotropic  laminate  (cantilever  simply  supported  plates). 
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Figure  121.  Weight  saving  for  different  cross  sections  whose  height  is  twice  the 
laminate  thickness  (one-dimensional  cantilever  plates). 


145 


T300/5208 


21  % 


WEIGHT  SAVING 
58  % 


66  % 


Figure  122.  Weight  saving  for  different  cross  sections  whose  height  is  ten  times  the 
laminate  thickness  (one-dimensional  simply  supported  plates). 


Figure  123.  Definition  of  coordinate  axis  for  two-dimensional  plates  subjected  to 
uniform  uniaxial  compression,  biaxial  compression,  shear  or  combined  load. 
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Figure  124.  Critical  load  and  optimum  angle  for  a  simply  supported  plate  subjected 
to  a  uniform  uniaxial  compression  load. 
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Figure  125.  Comparison  between  critical  loads  of  optimum  configuration  and 
quasi -isotropic  laminate  (simply  supported  plates  subjected  to  a  uniform  uniaxial 
compression  load). 
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Figure  126.  Comparison  of  critical  loads  for  cross-ply  and  quasi-isotropic  laminate 
for  a  square  plate  (simply  supported  plates  subjected  to  a  uniform  uniaxial 
compression  load). 
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Figure  127.  Comparison  between  weights  of  aluminium  and  the  optimum 
configuration  of  T300/N5208  (simply  supported  plates  subjected  to  a  uniform 
uniaxial  compression  load). 
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Figure  1 28.  Critical  load  and  optimum  angle  for  a  simply  supported  plate  subjected 
to  a  uniform  biaxial  compression  load  (Ny=Nx/2). 
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Figure  129.  Comparison  between  critical  loads  of  optimum  configuration  and 
quasi-isotropic  laminate  (simply  supported  plates  subjected  to  a  uniform  biaxial 
compression  load  (Ny=Nx/2)). 


nmmi 


Figure  130.  Comparison  of  critical  loads  for  cross-ply  and  quasi-isotropic  laminate 
for  a  square  plate  (simply  supported  plates  subjected  to  a  uniform  biaxial 
compression  load  (Ny=Nx/2)). 


nuiuu  20 

1.5 

timtttt  10 

RELATIVE  WEIGHT 
RESPECT  TO  0  5 

[ <t>  /-<t>  ]  opt 

o.o 

0  2  4  6  8  10  12 

_a_ 

b 

Figure  131.  Comparison  between  weights  of  aluminium  and  the  optimum 
configuration  of  T300/N5208  (simply  supported  plates  subjected  to  a  uniform 
biaxial  compression  load  (Ny=Nx/2)). 
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Figure  132.  Critical  load  and  optimum  angle  for  a  simply  supported  plate  subjected 
to  a  uniform  biaxial  compression  load  (Ny=Nx). 


Figure  13 3. Comparison  between  critical  loads  of  optimum  configuration  and  quasi- 
isotropic  laminate. (simply  supported  plates  subjected  to  a  uniform  biaxial 
compression  load  (Ny=Nx)). 
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Figure  134.  Comparison  of  critical  loads  for  cross-ply  and  quasi-isotropic  laminate 
for  a  square  plate  (simply  supported  plates  subjected  to  a  uniform  biaxial 
compression  load  (Ny=Nx)). 
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Figure  135.  Comparison  between  weights  of  aluminium  and  the  optimum 
configuration  of  T300/N5208  (simply  supported  plates  subjected  to  a  uniform 
biaxial  compression  load  (Ny=Nx)). 


0  2  4  6  8  10  12 

_0_ 

b 

Figure  1 36.  Critical  load  and  optimum  angle  for  a  simply  supported  plate  subjected 
to  a  uniform  biaxial  compression  load  (Nys=2*Nx). 
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Figure  137.  Comparison  between  critical  loads  of  optimum  configuration  and 
quasi-isotropic  laminate  (simply  supported  plates  subjected  to  a  uniform  biaxial 
compression  load  (Ny=2*Nx)). 
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Figure  138.  Comparison  of  critical  loads  for  cross-ply  and  quasi-isotropic  laminate 
for  a  square  plate  (simply  supported  plates  subjected  to  a  uniform  biaxial 
compression  load  (Ny=2*Nx)). 
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Figure  139.  Comparison  between  weights  of  aluminium  and  the  optimum 
configuration  of  T300/N5208  (simply  supported  plates  subjected  to  a  uniform 
biaxial  compression  load  (Ny=2*Nx)). 
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Figure  140.  Critical  load  and  optimum  angle  for  a  clamped  plate  subjected  to  a 
uniform  uniaxial  compression  load. 
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Figure  141.  Comparison  between  critical  loads  of  optimum  configuration  and 
quasi-isotropic  laminate  (clamped  plates  subjected  to  a  uniform  uniaxial 
compression  load). 
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Figure  143.  Comparison  between  weights  of  aluminium  and  the  optimum 
configuration  of  T300/N5208  (clamped  plates  subjected  to  a  uniform  uniaxial 
compression  load). 


0  2  4  6  8  10  12 

0 


b 

Figure  144.  Critical  load  and  optimum  angle  for  a  clamped  plate  subjected  to  a 
uniform  biaxial  compression  load  (Ny=Nx/2). 
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Figure  145.  Comparison  between  critical  loads  of  optimum  configuration  and 
quasi-isotropic  laminate  (clamped  plates  subjected  to  a  uniform  biaxial 
compression  load  (Ny=Nx/2)). 
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Figure  146.  Comparison  of  critical  loads  for  cross-ply  and  quasi-isotropic  laminate 
for  a  square  plate  (clamped  plates  subjected  to  a  uniform  biaxial  compression  load 
(Ny=Nx/2)). 
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Figure  147.  Comparison  between  weights  of  aluminium  and  the  optimum 
configuration  of  T300/N5208  (clamped  plates  subjected  to  a  uniform  biaxial 
compression  load  (Ny=Nx/2)). 
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Figure  148.  Critical  load  and  optimum  angle  for  a  clamped  plate  subjected  to  a 
uniform  biaxial  compression  load  (Ny=Nx). 
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Figure  149.  Comparison  between  critical  loads  of  optimum  configuration  and 
quasi-isotropic  laminate  (clamped  plates  subjected  to  a  uniform  biaxial 
compression  load  (Ny=Nx)). 
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Figure  150.  Comparison  of  critical  loads  for  cross-ply  and  quasi -isotropic  laminate 
for  a  square  plate  (clamped  plates  subjected  to  a  uniform  biaxial  compression  load 
(Ny=Nx)). 
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Figure  151.  Comparison  between  weights  of  aluminium  and  the  optimum 
configuration  of  T300/N5208  (clamped  plates  subjected  to  a  uniform  biaxial 
compression  load  (Ny=Nx)). 
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Figure  152.  Critical  load  and  optimum  angle  for  a  clamped  plate  subjected  to  a 


uniform  biaxial  compression  load  (Ny=2*Nx). 


Figure  153.  Comparison  between  critical  loads  of  optimum  configuration  and 
quasi-isotropic  laminate  (clamped  plates  subjected  to  a  uniform  biaxial  compression 


Figure  154.  Comparison  of  critical  loads  for  cross-ply  and  quasi-isotropic  laminate 
for  a  square  plate  (clamped  plates  subjected  to  a  uniform  biaxial  compression  load 
(Ny=2*Nx)). 
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Figure  155.  Comparison  between  weights  of  aluminium  and  the  optimum 
configuration  of  T300/N5208  (clamped  plates  subjected  to  a  uniform  biaxial 
compression  load  (Ny=2*Nx)). 
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Figure  156.  Critical  load  and  optimum  angle  for  a  simply  supported  plate  subjected 
to  a  uniform  shear  load. 


Figure  157.  Comparison  between  critical  loads  of  optimum  configuration  and 
quasi-isotropic  laminate  (simply  supported  plates  subjected  to  a  uniform  shear 


Figure  158.  Comparison  of  critical  loads  for  cross-ply  and  quasi-isotropic  laminate 
for  a  square  plate  (simply  supported  plates  subjected  to  a  uniform  shear  load). 


I 


b 


Figure  159.  Comparison  between  weights  of  aluminium  and  the  optimum 
configuration  of  T300/N5208  (simply  supported  plates  subjected  to  a  uniform 
shear  load). 
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Figure  160.  Critical  load  and  optimum  angle  for  a  clamped  plate  subjected  to  a 
uniform  shear  load. 
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Figure  161.  Comparison  between  critical  loads  of  optimum  configuration  and 
quasi-isotropic  laminate  (clamped  plates  subjected  to  a  uniform  shear  load). 
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Figure  162.  Comparison  of  critical  loads  for  cross-ply  and  quasi-isotropic  laminate 
for  a  square  plate  (clamped  plates  subjected  to  a  uniform  shear  load). 
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Figure  163.  Comparison  between  weights  of  aluminium  and  the  optimum 
configuration  of  T300/N5208  (clamped  plates  subjected  to  a  uniform  shear  load). 
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Figure  164.  Buckling  parameters  for  simply  supported  plates  subjected  to 
combined  uniform  uniaxial  compression  and  shear. 
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Figure  165.  Buckling  parameters  for  clamped  plates  subjected  to  combined 
uniform  uniaxial  compression  and  shear. 
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